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Retelle,  John  Powers,  Jr.  (Ph,D.,  Aerospace  Engineering  Sciences) 
Unsteady  Boundary  Layer  Flow  Reversal  in  a Longitudinally  Oscillating 
Flow 

Thesis  directed  by  Assistant  Professor  Donald  A.  Kennedy. 

This  experiment  was  conducted  to  provide  information  which 
would  lead  to  a better  understanding  of  the  onset  of  unsteady  flow 
reversal  in  attached  and  separated  boundary  layers.  It  is  shown 
that  boundary  layer  flow  reversal  on  a stationary  airfoil  in  an 
oscillating  free  stream  responds  to  changes  in  the  flow  unsteadiness 
in  a manner  similar  to  the  flow  reversal  noted  during  the  onset  of 
dynamic  stall  on  an  airfoil  oscillating  in  pitch.  Measurements  of 
boundary  layer  flow  reversal  made  in  a wind  tunnel  on  a stationary 
NACA  0012  airfoil  immersed  in  a longitudinally  oscillating  subsonic 
flow  are  described.  A split-film  sensor  provided  flow  reversal 
data  in  regions  above  the  airfoil  surface  throughout  the  period  of 
flow  oscillation.  At  angles  of  attack  below  the  steady-flow  stall 
angle  of  attack  and  for  small  amplitude  oscillations,  regions  of 
boundary  layer  fl^w  reversal  near  the  leading  edge  were  observed  to 
vary  nominally  in  phase  with  the  time-varying  external  velocity. 

Ho wever,  a large  increase  in  the  oscillation  amplitude  produced  a 
phase  lead  and  an  increased  downstream  extent  of  the  flow  reversal 
region.  Variations  in  the  oscillation  frequency  produced  only 
small  changes  in  the  shape  of  the  flow  reversal  region.  For  angles 
of  attack  greater  than  the  steady-flow  stall  angle,  the  growth  and 
decay  of  reversed  flow  regions  were  observed  to  be  stronger  functions 
of  the  velocity  oscillation  frequency  and  amplitude.  Increases  in 
frequency  were  observed  to  eliminate  flow  reversal  near  the  trailing 


i v 

edge.  Increases  in  oscillation  amplitude  produced  larger  reversal 
regions  which  nominally  lagged  the  phase  of  the  free  stream  velocity. 
Although  the  unsteady  wind  tunnel  was  adjusted  to  reduce  harmonic 
content,  the  presence  of  higher  harmonics  in  the  power  spectrum  of 
the  oscillating  free  stream  velocity  was  observed  to  produce  flow 
reversals  which  led  the  phase  of  the  free  stream  velocity.  The 
unsteady  effects  noted  above  are  attributed  to  both  the  temporal 
and  spatial  gradients  which  appear  in  the  unsteady  boundary  layer 
equations.  Examination  of  the  dimensionless  form  of  the  equations 
of  motion  reveals  that  variations  in  the  oscillation  frequency 
influences  only  the  temporal  gradients,  while  changes  in  the  oscilla- 
tion amplitude  can  produce  strongly  coupled  effects  on  both  the 
temporal  and  spatial  gradients.  Data  obtained  describing  the 
influence  of  the  unsteady  flow  parameters  on  flow  reversal  may  be 
applicable  to  the  initial  stages  of  dynamic  stall. 

This  abstract  is  approved  as  to  form  and  content.  I recommend  its 
publication. 

Signed'^  avv.uJ\Jl  1^*  - 

Faculty  member  in  charge  of  thesis 
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ABSTRACT 

Measurements  of  boundary  layer  flow  reversal  made  in  a wind  tunnel 
on  a stationary  NACA  0012  airfoil  immersed  in  a longitudinally  oscillating 
subsonic  flow  are  described.  A split-film  sensor  provided  flow  reversal 
data  in  regions  above  the  airfoil  surface  throughout  the  period  of  flow 
oscillation.  At  angles  of  attack  below  the  steady-flow  stall  angle  of 
attack  and  for  small  amplitude  oscillations,  regions  of  boundary  layer 
flow  reversal  near  the  leading  edge  were  observed  to  vary  nominally  in 
phase  with  the  time-varying  external  velocity.  However,  an  increase  in 
the  oscillation  amplitude  produced  a phase  lead  and  an  increased  down- 
stream extent  of  the  flow  reversal  region.  Variations  in  the  oscillation 
frequency  produced  only  small  changes  in  the  shape  of  the  flow  reversal 
region.  For  angles  of  attack  greater  than  the  steady-flow  stall  angle, 
the  growth  and  decay  of  reversed  flow  regions  were  observed  to  be  stronger 
functions  of  the  velocity  oscillation  frequency  and  amplitude.  Increases 
in  oscillation  amplitude  produced  larger  reversal  regions  which  nominally 
lagged  the  phase  of  the  free  stream  velocity.-  Although  the  unsteady  wind 
tunnel  was  adjusted  to  reduce  harmonic  content,  the  presence  of  higher 
harmonics  in  the  power  spectrum  of  the  oscillating  free  stream  velocity 
was  observed  to  produce  flow  reversals  which  led  the  phase  of  the  free 
stream  velocity.  It  is  shown  that  boundary  layer  flow  reversal  on  a 
stationary  airfoil  in  an  oscillating  free  stream  responds  to  changes  in 
the  flow  unsteadiness  in  a manner  similar  to  the  flow  reversal  observed 
during  the  onset  of  dynamic  stall  on  an  airfoil  oscillating  in  pitch. 
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1,  INTRODUCTION 


1,1  Statement  of  the  Problem 

This  experiment  was  conducted  to  provide  information  which 
would  lead  to  a better  understanding  of  the  onset  of  unsteady  flow 
reversal  in  attached  and  separated  boundary  layers.  Regions  of  flow 
reversal  are  produced  by  unsteady  flow  effects  and  are  often  encoun- 
tered in  flows  associated  with  maneuvering  aircraft  and  retreating 
helicopter  blades.  These  two  dynamic  applications  involve  rapid 
changes  of  both  angle  of  attack  and  flow  velocity.  Previous  experi- 
mental investigations  have  considered  these  parametric  variations 
separately.  Experiments  conducted  on  stationary  airfoils  in 
oscillatory  flows  which  are  described  below  have  involved  regions  of 
reversed  flow  imbedded  in  attached  boundary  layers.  Other  experi- 
ments concerned  with  dynamic  stall  phenomena  have  shown  that  airfoils 

I 

oscillating  in  pitch  may  produce  unsteady  values  of  lift  and  pitching 
moment  which  greatly  exceed  those  obtained  for  a stationary  airfoil. 

In  both  these  cases,  the  flow  unsteadiness  produced  an  unseparaced 
flow  reversal  region  on  the  airfoil  upper  surface  which  subsequently 
formed  a recirculating  flow  region  near  the  leading  edge.  The  only 
quantitative  measurements  of  the  ons:t  of  the  flow  reversal  region 
have  been  on  the  airfoil  surface  due  to  the  limitations  of  data 


collection  transducers. 
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^ 

Regions  of  unsteady  boundary  layer  flow  reversal  were 
measured  above  the  stationary  airfoil  which  were  related  to  those 
measured  at  the  surface  of  the  oscillating  airfoil.  However,  the 
spatial  extent  of  the  flow  reversal  above  the  surface  could  be 
measured  with  greater  ease  and  accuracy  above  the  fixed  airfoil. 
Also,  this  experiment  allowed  for  an  investigation  of  the  parameters 
of  the  unsteady  flow  field,  without  strong  coupling  with  the 
inertial  effects  of  the  airflow  about  the  oscillating  airfoil. 
Variations  in  the  frequency  and  amplitude  of  the  unsteady  velocity 
oscillations  were  observed  to  have  a strong  influence  on  the 
boundary  layer  flow  reversal.  This  influence  was  found  to  be 
similar  to  the  unsteady  effects  noted  with  the  dynamic  stall  of 
oscillating  airfoils. 

This  experiment  was  motivated  by  a desire  to  investigate  the 
feasibility  of  constructive  applications  of  dynamic  effects  for 
maneuvering  flight.  Aircraft  have  c’assically  been  designed  using 
aerodynamic  characteristics  derived  from  steady  flow  considerations 
since  experimental  data  were  not  available  to  describe  the  aero- 
dynamics of  unsteady  flow.  For  example,  the  solutions  to  the 
aircraft  equations  of  motion,  while  dynamic  in  themselves,  have  been 
generated  using  forcing  functions  based  on  steady  stability  deriva- 
tives. These  derivatives  have  been  calculated  using  wind  tunnel 
data  collected  at  fixed  values  of  velocity  and  angle  of  attack.  A 
more  thorough  understanding  of  the  unsteady  variations  of  lift  and 
pitching  moment  during  both  dynamic  stall  and  dynamic  velocity 
variation  might  allow  aircraft  to  be  designed  to  control  and  utilize 
the  high  lift  transients  generated  by  this  dynamic  process  for 


improved  maneuvering  capabilities.  Properly  configured  aircraft 
might  briefly  fly  outside  the  static  flight  envelope  which  has 
classically  been  used  to  define  aircraft  flight  limits. 
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1.2  Oscillating  Freestrearn  Velocity  Experiments 

The  boundary  layer  flow  reversals  which  often  initiate  the 
dynamic  stall  process  were  difficult  to  measure  in  moving-airfoil 
experiments.  Flow  sensors  could  not  easily  be  positioned  above  the 
moving  airfoil.  Flow  reversal  data  have  been  more  conveniently 
collected  on  a stationary  airfoil  immersed  in  a longitudinally 
oscillating  flow.  Such  experiments,  described  below,  have  proven 
more  suitable  for  describing  the  unsteady  effects  which  may  provide 
the  mechanisms  which  lead  to  dynamic  stall.  However,  the  global 
effects  of  the  classical  dynamic  stall,  namely  the  separation  of  a 
leading  edge  vortex,  and  the  large  hysteresis  effects  on  lift,  were 
not  seen  with  this  experimental  environment. 

A typical  reversed  flow  region  observed  above  the  fixed 
airfoil  at  one  phase  point  of  the  freestrearn  velocity  oscillation 
is  shown  in  Figure  1.  A thin  layer  of  reversed  flow  was  imbedded 
within  the  boundary  layer,  but  flow  remained  attached  to  the  airfoil 
surface.  This  unsteady  flow  result  was  markedly  different  from  the 
steady  flow  case  where  flow  reversal  (or  zero  skin  friction)  has 
defined  the  onset  of  boundary  layer  separation.  A typical  evolution 
of  the  reversed  flow  region  during  an  oscillation  period  is  shown 
in  Figure  2.  Flow  reversal  was  observed  to  start  near  the  leading 
edge  at  approximately  the  n-radian  phase  point  (Figure  2b)  in  the 
oscillation  cycle.  As  the  magnitude  of  the  free  stream  velocity 
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UNSTEADY  FLOW  REVERSAL  ON  STATIONARY  NACA  0012  AIRFOIL 


(a)  No  Flow  Reversal  Present 


(b)  Reversal  Region  Near 
Leading  Edge 


(c)  Reversal  Region  Extends 
to  Trailing  Edge 


(d)  Reversal  Region  Shrinks 
Upstream  Toward  Leading 
Edge 


(e)  No  Reversal  Region  Present 


Figure  2 

FLOW  REVERSAL  EVENTS  ON  STATIONARY  NACA  0012  AIRFOIL 
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decreased  (Figure  2c),  the  reversal  region  extended  further  down- 
stream with  a simultaneous  increase  in  height  near  the  leading 
edge.  Flow  reversal  was  generally  observed  to  be  absent  for  the 
first  portion  of  the  velocity  oscillation  period  (0  _<  wt  _<  tt). 

The  work  of  Saxena,  Fejer,  and  Morkovin^  describes  wind 
tunnel  measurements  of  boundary  layer  profiles  and  surface  pressure 
distributions  on  a stationary  NACA  0012  airfoil  in  an  oscillating 
flow.  Data  were  collected  using  hot-wire  anemometers,  surface- 
mounted  film  pressure  gages,  and  flow  visualization  using  tufts. 

Flow  reversals  were  not  quantitatively  measured  because  of  the 
limitations  of  the  single-element  hot-wires  used.  Their  experiments 
showed  the  boundary  layer  flow  to  be  quasi-steady,  that  is,  in 
phase  with  the  time-varying  external  velocity,  for  airfoil  angles  of 
attack  below  the  steady  flow  stall  angle.  For  larger  angles  of 
attack,  increases  in  the  free  stream  velocity  oscillation  frequency 
were  observed  to  cause  large  departures  from  quasi-steady  behavior. 
No  measurements  were  made  of  the  effect  of  changing  the  amplitude 
of  the  free  stream  velocity  oscillations. 

The  data  described  above  are  supplemented  and  amplified  by 
the  results  of  the  experiment  presented  below.  Boundary  layer  flow 
reversal  boundaries  were  measured  using  a split-film  sensor  in 
regions  at  and  above  the  upper  surface  of  a NACA  0012  airfoil.  Both 
the  free  stream  velocity  oscillation  amplitude  and  frequency  were 
varied  to  produce  unsteady  effects  on  the  boundary  layer  flow 
reversals.  The  departures  from  quasi-steady  behavior  were  found  to 
depend  on  oscillation  frequency  in  a similar  manner  to  that  noted 
by  Saxena,  Fejer,  and  MorkovinJ  However,  the  results  described 


here  are  based  upon  quantitative  measurements  of  flow  reversal 
rather  than  inference  from  surface  pressure  measurements.  These 
results  also  show  the  strong  influence  of  changes  in  the  velocity 
oscillation  magnitude.  The  effects  of  large  velocity  oscillation 
amplitudes  were  observed  in  both  the  departures  from  quasi-steady 
behavior  and  the  large  changes  in  the  shape  and  size  of  the  flow 
reversal  region. 

The  scheme  of  oscillating  the  longitudinal  velocity  about  a 

2 3 

mean  was  used  in  the  analytical  studies  of  Lighthill,  Lin,  and 
Nickerson.^  For  the  boundary  layer  above  a semi-infinite  flat 
plate,  the  analyses  indicated  that  both  the  wall  shear  and  boundary 
layer  thickness  were  functions  of  the  reduced  frequency  of  velocity 
oscillation.  The  reduced  frequency,  which  resembles  the  Strouhal 
number  in  form,  depends  on  the  oscillation  frequency  and  the  mean 
velocity.  For  low  frequencies,  the  boundary  layer  responded  to  the 
flow  unsteadiness  in  a quasi-steady  ■ anner,  i.e.,  in  phase  with  the 
free  stream  flow.  For  high  frequencies,  the  analysis  predicted 
that  the  unsteady  effects  were  confined  to  a thin  shear  layer  near 
the  wall  and  that  shear  response  phase  in  this  layer  led  the  free 
stream  by  45  degrees  for  laminar  flow  and  by  about  10-15  degrees 

5 

for  turbulent  flow.  Experimental  data  collected  by  Karlsson, 

6 7 8 

Hill  and  Stenning,  Cousteix,  et.  al . and  Patel  confirmed  these 
predictions  for  both  laminar  and  turbulent  boundary  layers  on  semi- 
infinite flat  plates.  Additional  confirmation  was  provided  by  the 

g 

numerical  studies  of  Patel  and  Nash  and  Tclionis,  Tsahalis,  and 
Werlc.^  Their  calculations  for  flat  plates  predicted  the  presence 
of  thin  regions  of  reversed  flow  ahead  of  the  point  of  separation 
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of  the  unsteady  boundary  layer.  The  flow  reversal  indicated  that, 
for  unsteady  flows,  there  was  no  correlation  between  zero  wall 
shear  and  separation.  The  effects  of  variations  of  the  unsteady 
flow  parameters  (velocity  oscillation  frequency  and  amplitude)  on 
the  flow  reversals  were  also  computed. 

Calculations  were  also  made  to  predict  the  point  of  separa- 
tion of  unsteady  boundary  layers  subjected  to  spatial  pressure 
gradients.  Numerical  results  were  used  by  Sears  and  Telionis^  to 
propose  the  unsteady  boundary  layer  separation  model  summarized  in 
Figure  3.  At  a particular  phase  point  in  the  velocity  oscillation 
cycle,  the  point  of  zero  surface  shear  is  followed  downstream  by  a 
thin  layer  of  reversed  flow.  The  model  shows  a thickened  reversed 
flow  region  further  downstream.  Sears  and  Telionis  contend  that 
the  start  of  this  aft  thickened  region  represents  the  separation 

point  of  the  unsteady  boundary  layer.  The  numerical  calculations 

12  13 

of  Nash  and  Scruggs,  and  Nash  also  predicted  an  increase  in  the 
displacement  boundary  layer  thickness  downstream  of  the  initial  flow 
reversal  region  but  made  no  predictions  of  flow  separation. 

Experiments  have  been  conducted  to  measure  the  influence  of 
oscillating  flows  on  both  boundary  layer  flow  reversals  and  flow 
separation.  The  influence  of  flow  unsteadiness  on  laminar  separa- 
tion from  a semi-infinite  flat  plate  was  measured  by  Despard  and 

1 4 

Miller.  The  point  of  separation  was  observed  to  occur  upstream 
of  the  steady  separation  point.  Changes  in  the  flow  Reynolds  number 
and  oscillation  frequency  greatly  influenced  the  separation  point 
location,  while  almost  no  effect  was  noted  due  to  changes  in  the 
velocity  oscillation  amplitude. 


Separation 


Figure  3 

FLOW  REVERSAL  AND  SEPARATION  MODEL 
(reference  14) 
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Other  experiments  measured  the  phase  shift  between  the 

unsteady  boundary  layer  response  and  the  oscillating  freestream. 

1 5 

Simpson  measured  flow  reversals  in  a turbulent  boundary  layer  on 
a flat  plate  with  a longitudinal  spatial  pressure  gradient.  Data 
were  collected  using  hot-wire  anemometers  and  a laser  flow  visuali- 
zation scheme.  The  phase  of  the  reversed  flow  region  was  noted  both 
upstream  and  downstream  of  the  flow  separation  point.  Upstream  of 
the  separation  point,  the  phase  of  the  reversed  flow  led  the  free 
stream  oscillations.  The  phase  lead  was  greatest  near  the  wall,  in 
the  viscous  sublayer,  and  decreased  with  greater  distance  from  the 
wall.  Downstream  of  the  separation  point,  the  phase  of  the  reversed 
flow  was  observed  to  lag  the  free  stream  by  a large  amount.  Simpson 
contended  that  the  lag  observed  in  the  separated  region  was  due  to 
the  larger  flow  eddies,  which  necessitated  a larger  time  for  the 
outer  shear  layer  to  provide  mass  to  the  inner  layer  of  flow 
reversal.  Also,  near  the  separation  point,  the  bursting  frequency 
of  the  turbulent  boundary  layer  was  observed  to  decrease.  As  this 
frequency  became  closer  to  the  free  stream  oscillation  frequency,  a 
greater  interaction  was  noted  between  the  boundary  layer  and  the 
free  stream. 

Kennison^  also  measured  flow  reversals  in  a turbulent 
boundary  layer  on  a flat  plate  with  an  adverse  spatial  pressure 
gradient.  However,  the  oscillating  free  stream  velocity  included  a 
traveling  wave  disturbance  which  convected  downstream  at  less  than 
the  mean  free  stream  velocity.  Hot-wire  anemometers  and  hot-film 
surface  gages  produced  data  which,  when  corrected,  measured  flow 
reversals.  For  this  type  of  flow  oscillation,  the  adverse  spatial 
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pressure  gradient  was  observed  to  produce  a phase  lag  in  the  flow 
reversal  region  which  increased  with  downstream  movement.  The 
greatest  magnitude  of  phase  lag  was  noted  in  the  outer  region  of 
the  boundary  layer.  The  thin  layer  of  flow  reversal  ahead  of  the 
flow  separation  point  had  little  effect  on  the  boundary  layer 
momentum  thickness,  indicating  that  there  was  little  vorticity 
transfer  into  the  free  stream.  The  phase  lag  of  the  flow  reversal 
was  probably  due  to  the  presence  of  the  traveling  wave  disturbance 
in  the  free  stream  velocity. 

1.3  Dynamic  Stall  Experiments 

Unsteady  flow  wind  tunnel  experiments  have  also  measured  the 
dynamic  stall  phenomena  caused  by  rapid  changes  in  airfoil  angle  of 
attack.  Previous  experiments  have  utilized  oscillating  airfoil 
motion  in  steady  wind  tunnel  flows.  McCroskey  and  Philippe17  used 
a NACA  0012  airfoil  oscillating  in  pitch  about  a fixed  angle  of 
attack  near  the  steady  stall  angle.  Surface  pressure  measurements 
were  made  to  determine  the  chanqr^  in  the  lift  and  pitching  moment 
as  the  angle  of  attack  was  vari  Results  showed  increased  values 
of  lift  and  pitching  moment  ovc  ie  steady  values,  as  well  as 
large  non-linear  hysterisis  effects  in  both  parameters.  Hot- film 
gages  and  hot-wire  anemometers  were  used  to  measure  flow  reversal. 
The  dynamic  stall  events  observed  on  a NACA  0012  airfoil  are 
presented  in  Figure  4 (reprinted  from  reference  19).  The  onset  of 
this  process  is  indicated  by  the  appearance  of  boundary  layer  flow 
reversal  in  Figure  4(b).  This  reversed  flow  is  observed  to  precede 
the  formation  of  a vortex-like  disturbance  near  the  leading  edge  in 


(a)  STATIC  STALL  ANGLE.  EXCEEDED 

(b)  FIRST  APPEARANCE  OF  FLOW 
REVERSAL  ON  SURFACE 


(c)  LARGE  EDDIES  APPEAR  IN 
BOUNDARY  LAYER 


(d)  FLOW  REVERSAL  SPREADS  OVER 


MUCH  OF  AIRFOIL  CHORD 


(e)  VORTEX  FORMS  NEAR 
LEADING  EDGE 


( f ) LIFT  SLOPE  INCREASES 


(g)  MOMENT  STALL  OCCURS 


(h)  LIFT  STALL  BEGINS 


( i ) MAXIMUM  NEGATIVE  MOMENT 
( i ) FULL  STALL 


(k)  BOUNDARY  LAYER  REATTACHES 
FRONT  TO  REAR 


( I ) RETURN  TO  UNSTALLED  VALUES 
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Figure  4(e).  The  subsequent  separation  of  this  vortex  from  the 
airfoil  surface  produced  a rapid  drop  in  the  pitching  moment 
(Figure  4(g)  - moment  stall)  and  then  a drop  in  normal  force 
(Figure  4(h)  - lift  stall).  Later  experiments  by  McCroskey,  et. 
al . **  and  Carr,  et.  alJ^  produced  data  showing  the  changes  in  the 
lift,  pitching  moment,  and  flow  reversals  produced  by  varying  the 
flow  Reynolds  number,  mean  angle  of  attack,  airfoil  leading  edge 
geometry,  and  the  airfoil  oscillation  frequency  and  amplitude. 

Their  flow  reversal  measurements  made  at  the  airfoil  surface  are  of 
interest  in  this  discussion.  A typical  result  is  presented  in 
Figure  5.  The  horizontal  axis  represents  nondimensional  chord, 
while  the  vertical  axis  represents  the  phase  and  magnitude  of  the 
sinusoidal  variation  of  angle  of  attack.  The  symbols  located  near 
a phase  of  tt/4  represent  the  start  of  reversed  flow  prior  to  flow 
separation.  The  symbols  near  a phase  of  5n/4  represent  the 
reattachment  of  the  boundary  layer  tr.  the  airfoil  surface.  These 
investigators  also  reported  that  boundary  layer  separation  occurred 
initially  near  the  trailing  edge  for  the  unmodified  NACA  0012 
airfoil.  A modification  to  decrease  the  radius  of  the  airfoil 
leading  edge  caused  separation  to  initially  occur  much  closer  to 
the  leading  edge. 

Airfoil  motion  was  produced  in  a different  manner  with  the 

experimental  apparatus  described  by  Rebont,  Maresca,  Favier  and 
20 

Valensi.  They  oscillated  the  airfoil  fore  and  aft  in  translation 
at  a fixed  angle  of  attack  in  steady  flow.  For  angles  of  attack 
above  the  steady  flow  stall  angle,  their  results  showed  strong 
unsteady  effects  on  the  skin  friction.  The  chronology  of  the 
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dynamic  stall  events  observed  was  similar  to  the  results  obtained 
for  the  airfoil  oscillating  in  pitch. 

1.4  Optimization  of  Oscillating  Free  Stream 

The  presence  of  higher  harmonic  content  in  the  power  spectrum 

of  the  free  stream  velocity  was  observed  to  produce  large  unsteady 

effects  in  the  response  of  the  flow  reversal  region.  In  order  to 

reduce  these  effects,  a preliminary  experiment  was  conducted 

concerning  the  harmonic  optimization  of  the  unsteady  velocity 

oscillations.  The  wind  tunnel  modification  described  below  was 

specifically  designed  to  produce  a sinusoidal  perturbation  to  the 

mean  velocity  which  was  free  from  unwanted  higher  harmonic  effects. 

A set  of  rotating  vanes  at  the  downstream  end  of  the  test  section 

was  designed  to  periodically  block  the  flow  to  the  axial  flow  fan 

located  further  downstream.  A nearly  constant  air  flow  was  provided 

to  the  fan  by  a bypass  mechanism  which  allowed  the  fan  to  operate 

at  a nearly  constant  pressure  rise,  improving  the  harmonic 

suppression  in  the  flow.  It  was  found  that  this  harmonic 

suppression  could  be  optimized  for  the  various  combinations  of 

velocity  oscillation  amplitude  and  frequency  used  in  this  experiment. 

Several  experiments  in  wind  tunnels  at  ether  institutions 

have  used  the  rotating  vane  principle  to  perturb  the  flow.  Two 

similar  tunnels  made  no  attempt  to  improve  the  suppression  of  higher 

harmonics.  An  apparatus  on  an  open-circuit  wind  tunnel  described 
21 

by  Despard  used  four  rotating  vanes  downstream  of  the  test  section 


to  produce  velocity  variations  of  up  to  40  percent  at  oscillation 
frequencies  of  up  to  240  Hz.  The  second  and  third  harmonics 
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contributed  23  percent  and  4 percent  of  the  velocity  power  spectrum, 

respectively.  An  almost  identical  apparatus  was  described  by 

22 

Morkovin,  Loehrke,  and  Fejer  which  perturbed  the  flow  in  a closed- 
circuit  wind  tunnel.  Resonances  were  developed  corresponding  to  the 
tunnel  circuit  length  and  also  to  the  distance  between  the  rotating 
vanes  and  the  tunnel  contraction  and  honeycomb.  Large  departures 
from  sinusoidal i ty  were  noted  for  large  velocity  oscillation  ampli- 
tudes. Recent  work  in  this  tunnel  by  Saxena,  Fejer,  and  Morkovin^ 
showed  plots  of  the  velocity  and  pressure  fluctuations  which 
indicate  a higher  harmonic  content  in  the  pressure  fluctuation  than 
in  the  longitudinal  velocity  variation. 

Several  different  schemes  have  been  tried  to  eliminate  higher 
harmonic  content  in  the  velocity  oscillations  of  these  rotating-vane 
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tunnels,  Karlsson  used  four  rotating  vanes  downstream  of  the  test 

section  of  an  open-circuit  wind  tunnel  to  produce  root-mean-square 

velocity  variations  of  34  percent  of  vane  rotation  frequencies  of 

up  to  48  Hz.  At  high  frequencies  the  higher  harmonics  were 

responsible  for  about  15  percent  of  the  total  root-mean-square 

amplitude  of  the  fluctuation.  In  an  effort  to  reduce  these 

harmonics,  the  lowest  vane  was  turned  at  a higher  rate,  but  no 

1 5 

mention  is  made  of  the  success  of  this  device.  Simpson  described 
an  apparatus  of  this  type  which  used  five  rotating  vanes  downstream 
of  the  test  section  of  an  open-circuit  tunnel.  A programmable  motor 
was  used  to  adjust  the  angular  velocity  of  the  rotating  vanes, 
thereby  optimizing  the  higher  harmonic  suppression.  The  resulting 
velocity  waveforms  appeared  sinusoidal,  but  no  harmonic  analysis 
was  presented. 
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An  interesting  variation  to  the  rotating  vane  device  was 

23 

described  by  Charnay  and  Mathieu.  Their  eight  rotating  vanes 
were  positioned  in  the  settling  chamber  ahead  of  the  screens  and 
honeycomb  sections  and  upstream  of  the  test  section  of  the  open- 
circuit  wind  tunnel.  They  found  that  the  smallest  harmonic  content 
occurred  when  the  vane  vortex-sheading  frequencies  matched  the 
tunnel  resonant  frequencies.  Turbulent  velocities  were  measured  to 
evaluate  the  effect  of  the  upstream-mounted  vanes.  It  was  observed 
that  optimum  sinusoidal ity  was  indicated  by  maxima  of  the  magnitudes 
of  the  longitudinal  and  vertical  turbulent  fluctuating  velocities, 
Uj1  and  u2',  and  also  by  minima  of  the  skewness  and  flatness  factors 
of  Uj'.  The  turbulence  intensity  in  the  test  section  downstream  of 
the  vanes  was  measured  to  be  on  the  order  of  0.7  percent. 

The  rotating  vane  apparatus  described  below  was  capable  of 

higher  harmonic  suppression  to  a greater  degree  than  the  devices 

described  above,  with  the  possible  exception  of  the  device 

1 5 

described  by  Simpson.  Careful  adjustment  of  the  bypass  massflow 
provided  regulation  of  the  higher  harmonics  in  the  test  section  for 
a given  tunnel  mean  velocity  and  vane  rotation  frequency.  Depending 
on  *he  vane  size  and  flow  conditions,  the  second  harmonic  power 
amplitude  could  be  suppressed  as  much  as  50  db. 


2,  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

2,1  Unsteady  Wind  Tunnel  Configuration 

This  experiment  was  conducted  in  the  two  foot  by  two  foot 
(0.61  meter  by  0.61  meter)  cross-section  low  speed  wind  tunnel  at 
the  University  of  Colorado,  It  was  an  open  return  suction-type 
wind  tunnel  with  the  addition  of  the  unsteady  flow  apparatus.  The 
unmodified  steady  flow  facility  was  described  in  detail  by 
Francis,^4  The  wind  tunnel  schematic  is  shown  in  Figure  6. 

Upstream  of  the  test  section  were  located  a plenum  with  a hexagonal 
honeycomb,  four  screens,  and  a sixteen-to-one  area  contraction.  The 
test  section  was  14  feet  (4,27  meters)  in  length  with  a 10  foot 
(3,05  meter)  long  plexiglass  wall  of  1/2  inch  thickness  and  is  shown 
in  Figures  7 and  8,  Downstream  of  the  test  section  were  the 
unsteady  flow  device,  a transition  section  to  the  fan,  the  fan,  and 
the  conical  diffuser. 

A periodic  disturbance  to  the  mean  flow  was  produced  by  a 
set  of  six  rotating  vanes  located  just  downstream  of  the  test 
section.  Figure  9 shows  these  vanes  from  upstream  in  a position  of 
partial  flow  resistance.  Three  sets  of  vanes  of  varying  width  were 
used  dixing  this  experiment  to  change  the  amount  of  maximum  tunnel 
blockage  and  thus  the  velocity  amplitude  variation,  N.  The  vanes 
are  described  in  Table  1 below  and  also  in  Figure  10.  They  were 
made  in  a modified  diamond  shape,  both  for  ease  of  construction  and 
for  minimal  pressure  and  friction  losses. 


Figure  8 

TEST  SECTION  SCHEMATIC 


TABLE  I 


ROTATING  VANE  DESCRIPTION 


Vane  Set 

Width 

(inch) 

Thickness 
( inch) 

Maximum 

Blockage 

Ratio 

Minimum 

Blockage 

Ratio 

Set  A 
(smal 1 ) 

2,000 

0,750 

0,500 

0,188 

Set  B 
(medium) 

3.495 

0.750 

0.874 

0.188 

Set  C 
(large) 

3.980 

0.750 

0.995 

0.188 

The  rotating  vanes  were  driven  by  a gear  train  powered  by  a General 
Electric  Stratitrol  II  3/4  horsepower  full -wave  direct-current  motor 
controlled  by  a Variac-equipped  direct-current  power  supply.  An 
8-to-l  gear  reduction  allowed  the  motor  to  run  smoothly  near  the 
middle  of  its  0-1759  revolution-per-minute  speed  range  which 
prevented  surging  at  low  speeds,  The  gear  drive  was  designed  to 
allow  the  vanes  to  counter-rotate  and,  thus,  produced  no  net  local 
lift  perturbation  to  the  flow, 

Figure  11  shows  a top  view  of  the  unsteady  flow  generator. 
This  view  shows  the  bypass  mechan'sm  which  allowed  air  to  enter  the 
fan  section  when  the  test  section  was  blocked  by  the  front,  or 
free  stream-controll ing,  vanes.  Three  rotating  vanes  on  the  top  of 
the  device  and  three  on  the  bottom  were  driven  by  the  same  gear 
train  which  also  drove  the  front  vanes.  The  vanes  on  the  top  and 
bottom  allowed  air  from  the  laboratory  to  enter  the  fan  section  at 
a varying  rate  which  balanced  the  disturbance  in  the  test  section. 
The  adjustment  of  the  bypass  air  flow  into  the  fan  section  was 


Figure  11 

BYPASS  VANES  AND  ADJUSTABLE  SLIDES 
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controlled  by  moveable  slides  which  are  shown  in  Figure  11,  These 
slides  were  set  to  predetermined  bypass  openings  based  on  the  flow 
conditions  desired  in  the  test  section.  Although  the  front  vanes 
were  changed  to  provide  changes  in  the  velocity  oscillation  ampli- 
tude, the  bypass  vanes  remained  at  a 2-inch  (5.08  cm)  width  (vane 
set  A),  This  arrangement  greatly  facilitated  the  changing  of  vane 
sets, 

A triggering  mechanism  based  on  rotating  vane  position  used 
a magnetic  pickup  mounted  on  a vane  shaft  to  produce  a voltage  pulse 
every  180  degrees  of  vane  rotation.  The  time  interval  between 
pulses  was  measured  by  a Computer  Measurements  model  605A  electronic 
counter.  This  enabled  the  half-period  of  vane  rotation  (or  total 
period  of  velocity  oscillation)  to  be  monitored  and  regulated  by 
adjustments  of  the  vane  rotation  speed.  Monitoring  of  the  velocity 
oscillation  period  showed  an  accuracy  of  + 0.005  seconds,  which 
represented  a frequency  accuracy  of  1.4  percent  at  a rotation 
frequency  of  1.43  Hz.  Another  circuit,  consisting  of  a regulated 
voltage  applied  to  a continuous-turn  potentiometer  mounted  on  a 
vane  shaft,  provided  a linear  output  voltage  corresponding  to  vane 
angle.  This  output  was  used  to  align  the  front  vanes  parallel  to 
the  flow  when  steady  conditions  were  desired  in  the  test  section. 

2.2  Airfoil  Configuration 

The  airfoil  used  for  this  ex,  eriment.  was  a NACA  0012  with  a 
10  inch  (25.4  cm)  chord,  A span  of  23,900  inches  (60.706  cm)  left 
a gap  of  0,050  inches  (0.127  cm)  between  the  wing  and  each  side  wall. 
The  solid  wing  was  milled  from  7075-T7  aluminum  to  a + 0.002  inch 
(0,005  cm)  tolerance  over  the  entire  profile.  The  wing  was  supported 
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at  the  quarter  chord  on  pins  which  protruded  through  oversized  holes 
in  the  wind  tunnel  walls.  The  pins  were  supported  from  outside  the 
tunnel  on  ball-bearings  by  a U-shaped  frame  placed  across  the 
I-beams  which  formed  the  top  of  the  test  section.  This  support 
served  to  isolate  the  airfoil  from  deflections  of  the  test  section 
side  walls  caused  by  the  unsteady  flow. 

The  airfoil  angle  of  attack  was  measured  by  aligning  a line 
scribed  along  the  airfoil  centerline  with  a leveled  protractor 
attached  to  the  plexiglass  test  section  wall.  Parallax  errors  were 
minimized  by  always  sighting  the  angle  of  attack  from  the  airfoil 
trailing  edge,  giving  values  of  angle  of  attack  accurate  to  within 
+ 0.2  degrees. 

2.3  Velocity  Measurements 

Free  stream  velocity  measurements  were  made  using  a constant- 
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temperature  hot-wire  anemometer  circuit  designed  by  Bank  and  shown 

in  Figure  12.  The  anemometer  output  voltage,  E,  can  be  related  to 
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the  free  stream  velocity  using  King's  law: 

E2  - E 2 = BU  n, 

o °° 

where  E is  a constant  voltage  level. 

0 3 

An  analog  circuit,  shown  in  Figure  13,  utilized  Motorola 

MCI  594  multipliers  to  linearize  the  anemometer  output.  Data 
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collected  by  Francis  in  this  facility  validated  the  use  of  this 


apparatus  by  showing  that  the  exponent  "n"  was  very  close  to  0.5 
and  that  the  coefficient  B was  a constant  for  flow  experiments  of 
this  type  using  air  and  hot-wire  probes  of  this  type.  It  was  also 


Figure  12 

HOT-WIRE  ANEMOMETER 
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shown,  using  a square  wave  test,  that  the  anemometer  frequency 
response  (+  3 db)  was  at  least  17  KHz.  This  response  was  recon- 
firmed for  the  hot  wires  used  in  the  present  experiment. 

The  hot-wire  probes  used  were  made  with  silver  coated 
Wollaston  wire  of  platinum  and  10  percent  rhodium  composition. 

The  wire  core  was  0.0001  inches  (.000254  cm)  in  diameter.  Electro- 
chemical etching  of  the  silver  coating  provided  a probe  sensing 
element  with  a ratio  of  length/diameter  of  approximately  300.  Both 
single-wire  and  x-wire  probes  were  used  for  data  collection. 
Linearized  velocity  output  from  the  x-wire  probe  was  resolved  into 

components  using  analog  addition  and  difference  circuits  shown  in 
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Figure  14.  A circuit  constructed  by  McMichael  , shown  in  Figure 
15,  was  used  for  determining  the  turbulent  (a.c.)  components  of  the 
x-wire  signals. 

A reference  for  fluid  velocity  measurement  was  provided  by  a 
pitot  static  tube,  located  at  the  downstream  end  of  the  test 
section.  It  was  connected  to  a Statham  model  PM5TC  differential 
pressure  transducer  with  a range  of  + 0.05  psid.  The  accurate 
measurement  of  the  steady  velocity  by  the  pressure  transducer  was 
essential  to  the  proper  calibration  of  the  hot-wire  linearizer 


output. 

The  free  stream  steady  turbulence  level, 


for  the  basic  wind  tunnel  was  measured  by  Francis* 


nominal 


value  of  0.06  percent  for  a free  stream  velocity  of  60  feet  per 
second  (18.29  meters  per  second).  Turbulence  intensity  measurements 
with  the  rotating  vane  device  installed  we re  made  using  a Bruel 


Figure  15 

C.  SIGNAL  CONDITIONER 
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and  Kjaer  type  2417  Random  Noise  Voltmeter  which  had  a bandwidth  of 


2-20,000  Hz.  With  the  large  front  vanes  (set  C)  fixed  at  the 

minimum  resistance  (full  open)  position,  the  steady  flow  turbulence 

level  was  measured  to  be  0.08  percent  at  a free  stream  velocity  of 

60  feet  per  second  (18.29  meters  per  second).  The  vanes  were  then 

rotated  at  3.33  Hz  and  the  turbulence  level  was  measured  to  be 

0.54  percent  at  a mean  flow  velocity  of  12  feet  per  second  (3.66 

meters  per  second).  To  determine  the  effect  of  unsteadiness  on  the 

higher  frequency  turbulence,  this  velocity  data  signal  was  filtered 

through  a 4 Hz  high-pass  filter.  The  higher  frequency  root-mean- 

square  turbulence  level  was  then  measured  to  be  0.05  percent. 

Reversals  in  flow  direction  of  the  longitudinal  velocity 

component  in  the  airfoil  boundary  layer  were  sensed  by  a model  1287V 
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split  film  sensor  made  by  TSI,  Inc.  The  sensor  geometry  is  shown 
in  Figure  16.  In  this  design,  two  electrically  independent  platinum 
films  were  plated  on  a single  quartz  fiber.  The  plane  of  the  split 
between  the  two  films  was  oriented  perpendicular  to  the  free  stream 
direction.  Each  film  segment  was  operated  independently  with  a 
separate  hot-wire  anemometer  circuit.  The  data  in  Figure  17  show 
the  square  of  the  average  of  the  two  anemometer  outputs  as  a 
function  of  the  square  root  of  velocity  with  the  plane  of  the  split 
between  the  probes  aligned  parallel  to  the  steady  flow.  The 
linearity  of  these  data  show  that  a relationship  similar  to  King's 
law  exists  between  anemometer  output  voltage  and  velocity.  The 
exact  formulation  of  this  relation  may  be  somewhat  altered  from 
King's  law,  but  the  linear  approximation  shown  in  Figure  17  is 
sufficient  over  the  velocity  range  measured  in  this  experiment. 

L - .. _ j 
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Measurements  by  TSI,  Inc.  showed  this  probe  to  be  very  insensitive 
to  errors  caused  by  both  fluid  velocities  along  the  axis  of  the 
quartz  fiber  and  by  the  effects  of  the  small  supporting  structures. 
The  frequency  response  of  this  probe  used  with  the  anemometer 
circuit  was  measured  using  the  square  wave  test  to  be  at  least  200 
Hz. 

A precision  three-dimensional  traversing  mechanism  allowed 
the  spatial  positioning  resolution  of  the  split  film  probe  to  be 
determined  to  within  + 0.005  inches  (0.013  cm)  in  horizontal  and 
axial  movement  and  to  within  +0.002  inches  (.005  cm)  in  vertical 
movement.  Motion  in  the  three  directions  was  provided  by  ball- 
bearing screw  drives  along  the  vertical  and  axial  axes  and  a rack- 
and-pinion  drive  along  the  horizontal  axis.  These  drives  were 
powered  by  miniature  direct  current  motors.  Distances  were  deter- 
mined with  geared  potentiometers  connected  to  the  drives  and  the 
amplifier  circuitry  shown  in  Figure  18.  Each  circuit  contained 
potentiometers  labeled  "span"  and  "reference  zero"  which  were  set  to 
provide  output  voltages  corresponding  to  position  at  each  end  of 
the  maximum  travel  of  the  probe.  The  potentiometer  marked  "position 
zero"  could  be  used  to  set  an  output  voltage  of  zero  at  any  inter- 
mediate position  of  the  probe  along  its  travel.  Probe  positions 
relative  to  the  airfoil  in  the  axial  and  vertical  axes  were  deter- 
mined using  the  airfoil  trailing  edoe  as  a baseline  reference 
position.  The  probe  position  at  the  trailing  edge  was  measured 
optically  after  each  change  of  angle  of  attack  using  a Gaertner 
optical  cathotometer . The  distances  of  the  probe  vertically  above 
the  airfoil  surface  and  horizontally  along  the  chord  during  data 
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collection  were  then  computed  using  programs  for  a Texas  Instrument 
TI-58  programmable  calculator  which  are  shown  in  Appendix  A. 

2.4  Flow  Reversal  Measurement 

Each  film  of  the  split-film  probe  was  connected  to  a separate 
hot-wire  anemometer  circuit.  For  preliminary  adjustment,  the  probe 
was  mounted  in  the  tunnel  with  the  plane  of  the  split  between  the 
films  parallel  to  the  flow.  With  a steady  free  stream  tunnel 
velocity  of  about  60  ft/sec  (18.29  meters  per  second),  the  ratio  of 
the  output  of  the  two  anemometer  circuits  was  adjusted  to  be 
approximately  one.  The  free  stream  velocity  was  then  changed  and 
the  ratio  of  the  output  voltages  was  readjusted  so  that  the  ratio 
remained  constant  with  changes  in  velocity.  Measurements  indicated 
that  this  ratio  remained  between  0.990  and  1.010  for  the  velocities 
encountered  in  this  experiment. 

The  split-film  probe  was  calibrated  for  flow  reversal  indica- 
tion using  a support  jig  mounted  on  the  plexiglass  test  section 
wall.  Using  this  device,  the  angle  between  the  probe  and  the  free 
stream  velocity  could  be  varied  from  outside  the  tunnel  and 
measured  geometrically  with  a protractor.  The  probe  was  turned  so 
that  the  plane  of  the  split  was  perpendicular  to  the  flow.  The 
output  voltages  from  the  two  anemometer  circuits  were  subtracted 
using  the  analog  difference  circuit  shown  in  Figure  14.  The  output 
of  this  circuit  was  then  measured  as  the  steady  velocity  was  varied 
from  about  60  feet  per  second  (18.3  meters  per  second)  to  zero  feet 
per  second.  This  procedure  was  then  repeated  with  the  probe  turned 
180  degrees.  A sample  calibration  curve  of  the  difference  circuit 


voltage  as  a function  of  steady  velocity  is  shown  in  Figure  19. 

Flow  reversal  was  attributed  to  a voltage  level  of  0.49  volts  for 
this  calibration.  The  probe  was  then  mounted  on  the  traverse 
mechanism  and  was  calibrated  for  forward  velocities  only.  The  data 
shown  in  Figure  19  show  that  the  threshold  voltage  indicating  flow 
reversal  remained  essentially  the  same  when  the  probe  was  moved  from 
the  jig  to  the  traverse.  This  threshold  voltage  v/as  not  zero,  since 
the  ratio  of  the  anemometer  output  voltages  set  above,  with  the 
plane  of  the  splits  parallel  to  the  flew,  was  never  exactly  equal 
to  one.  This  threshold  reversal  voltage  was  rechecked  prior  to 
data  collection  and  showed  excellent  repeatability.  The  total  drift 
of  the  reversal  threshold  voltage  was  observed  to  be  + .05  volts 
over  seven  consecutive  days  of  data  collection.  The  probe  response 
shown  is  non-linear  since  the  linearizer  circuits  were  not  used. 

This  was  done  to  improve  accuracy,  since  the  low  velocities  near 
flow  reversal  would  produce  linearizer  voltages  within  the  noise 
levels  of  the  multiplier  elements. 

The  split  film  probe  was  used  to  indicate  flow  reversal,  that 
is,  the  point  when  the  sign  of  the  longitudinal  velocity  becomes 
negative.  The  actual  magnitude  of  the  reversed  velocity  v/as  not 
measured  because  of  the  inability  of  the  split  film  probe  to  resolve 
the  reversed  velocity  components  accurately.  Characteristics 
provided  by  the  probe  manufacturer  showed  that  flows  aligned  closely 
with  the  plane  of  the  split  could  be  accurately  resolved  into 
components.  However,  the  measured  magnitudes  of  the  components  of 
flows  perpendicular  to  the  plane  of  the  split,  such  as  the  reversed 
flows,  were  shown  to  be  in  error  by  up  to  10  percent.  For  this 
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reason,  only  the  onset  of  flow  reversal  was  measured.  The  error 
band  on  the  probe  indication  of  flow  reversal  onset  depended 
strongly  on  the  ensemble  averaging  technique  used  in  this  experiment 
and  will  be  discussed  below. 

2.5  Data  Processing 

The  complete  signal  processing  apparatus  is  shown  in  Figure 

20.  The  heart  of  the  data  analysis  was  a Hewlett-Packard  model 
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5452B  Fourier  Analyzer  System.  This  device  is  a digital  instru- 
ment capable  of  analysis  of  analog  input  data  in  either  the  time 
domain  or  the  frequency  domain.  The  analog  input  amplitude  range 
extended  to  10  volts  with  a maximum  frequency  response  of  25  KHz. 

The  system  included  an  anal og-to-digi tal  converter,  a digital 
computer  with  an  8192-word  memory,  and  a digital  oscilloscope.  A 
hardwired  function  keyboard  was  provided  to  initiate  analyzer 
operations.  It  v/as  also  used  for  programmable  scaling  of  the 
digitally  converted  input  and  the  programming  of  repetitive  opera- 
tions . 

The  analog  input  data  was  digitized  so  that  each  word 
represented  a value  of  the  input  at  a given  point  in  time  for  a 
time  domain  input,  or  a value  of  a spectrum  at  a given  frequency 
for  frequency  domain  input.  Each  data  word  was  a 16-bit  binary 
number  representing  an  integer  from  -32768  to  +32767.  Groups  of 
data  words  were  collected  in  data  blocks  containing  a quantity  of 
data  words  equal  to  any  power  of  2 from  2^  = 128  to  2^  = 1024. 
Computer  operations  were  then  carried  out  on  the  data  blocks,  in 
either  the  frequency  or  time  domain,  as  appropriate.  The  sizing  of 
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the  data  blocks  in  both  domains  represented  a compromise  between  the 
accuracy  required  and  the  period  of  time  required  to  analyze  all  the 
data  generated  by  this  experiment.  In  the  time  domain,  the  maximum 
size  of  the  time  sample,  T,  was  chosen  to  be  1.0  seconds  and  the 
number  of  data  words  in  the  data  block.  N,  was  selected  to  be  256. 
This  choice  allowed  the  analyzer  output  to  be  generated  quickly 
using  the  relationship 

T = N[At] . 

This  provided  an  interval  of  0.0039  seconds  between  time  domain  data 
words  which  was  found  to  give  sufficient  accuracy.  In  the  frequency 
domain,  the  maximum  frequency  was  set  at  10  Hz  after  an  inspection 
of  higher  frequencies  revealed  no  power  spectral  data  of  signifi- 
cance to  this  experiment  (turbulence  neglected).  Using  the  same 
block  size  of  256  words,  the  relationship 


showed  that  the  frequency  resolution  was  0.078  Hz  between  frequency 
domain  data  words.  The  number  of  frequency  data  words  per  block  was 
N/2,  or  half  the  block  size,  since  the  frequency  information  was 
broken  into  real  and  imaginary  displays. 

An  important  feature  of  the  Fourier  Analyzer  was  its  ability 
to  perform  ensemble  averages  using  either  time  or  frequency  domain 
data.  The  concept  of  the  ensemble  average  may  be  understood  by 
considering,  for  example,  velocity  measurements  taken  in  the  time 
domain.  The  velocity  ensemble  average  was  determined  by  averaging 
over  a collection  of  data  blocks,  where  each  block  had  a common 


44 


phase  with  respect  to  the  varying  free  stream  velocity.  The 
velocity  field  U^(x,t;  p),  as  a function  of  a position  vector  x and 
a time  t,  was  given  an  ensemble  counter  parameter  p.  The  ensemble 


average  was  then: 


<Ui (x,t;  p)>  = 1 im 

m-*» 


,<(=,t|^. 
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Writing  the  velocity  as  a mean  plus  a deviation  about  a mean: 

U.j(x,t)  =<Ui(x,t)>+  u.'(x.t) 

it  can  be  seen  that  the  process  of  ensemble  averaging  eliminates  the 
deviations  about  the  mean,  since  u /r  = 0.  For  this  experiment,  the 
number  of  data  blocks,  m,  used  for  each  ensemble  average  was  chosen 
to  be  large  enough  so  that  the  deviations  from  the  mean  tended  to 
zero,  but  small  enough  so  that  the  period  of  time  required  for  data 
collection  was  not  excessive.  Experience  showed  that  ensemble 
averages  of  fifty  blocks  for  the  flow  field  analysis  were  sufficient 
to  separate  the  mean  velocity  from  the  input  velocity  data.  When 
the  flow  reverse1  data  was  subsequently  collected,  only  twenty  data 
blocks  were  averaged.  The  use  of  this  smaller  number  of  blocks 
reduced  the  total  time  required  tj  calculate  the  many  averages 
necessary  to  describe  the  flow.  This  smaller  number  was  found  to 
be  sufficient  to  eliminate  the  deviations  from  the  mean.  The  phase 
of  the  data  blocks  used  in  averaging  was  held  constant  by  the 
trigger  mechanism  connected  to  the  rotating  vanes.  The  trigger 


voltage  pulse  was  used  to  start  the  input  of  data  into  each  block. 
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2.6  Sources  of  Error 

The  error  band  on  the  flow  reversal  indication  measured  by 
the  split  film  probe  depended  on  the  ensemble  averaging  process. 

The  reversal  threshold  voltage  indicated  in  Figure  19  is  actually 
an  average  value  based  on  an  ensemble  of  20  data  blocks.  The  root- 
mean-square  deviation  of  the  threshold  voltages  from  each  block, 

V.;  from  the  ensemble  averaged  threshold  voltage  of  <V>  = 0.49 
volts;  was  determined  using  the  expression 
20 

y]  [(v.  - < v >)2i  v v 

L'  ^ ' ] ] average 

For  the  data  used  in  Figure  19,  this  root-mean-square  error  ratio 
was  calculated  to  be  1.1  percent  for  an  average  threshold  voltage 
of  0.49  volts.  Through  the  calibration  shown  in  Figure  18,  this 
corresponds  to  an  error  in  the  reversal  velocity  of  less  than  one- 
half  foot  per  second  (0.15  meters  per  second). 

In  addition  to  the  measurement  errors  already  mentioned, 
there  was  some  concern  over  tunnel -induced  three-dimensional 
effects.  At  high  oscillation  frequencies,  the  time-varying  velocity 
due  to  the  large  vanes  (set  C)  caused  considerable  physical  deflec- 
tion of  the  plexiglass  test  section  wall.  Unsteady  velocity 
measurements,  made  at  several  spanwise  locations  in  the  wing 
boundary  layer,  showed  that  these  three-dimensional  effects  were 
negligible.  The  measured  vertical  and  longitudinal  velocity 
components  varied  by  less  than  + 0.2  percent  along  the  wing  span 
outside  of  the  side-wall  boundary  layer. 
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Ambient  temperature  variations  were  minimized  by  operating 
the  wind  tunnel  at  or  below  about  60  feet  per  second  (18.29  meters 
per  second)  for  at  least  two  hours  before  starting  data  collection. 
This  allowed  the  temperature  of  the  laboratory  environment  to 
stabilize  at  approximately  70  degrees  Fahrenheit  (21  degrees 
Celsius).  Errors  due  to  electronic  drift  were  then  negligible. 
Calibrations  of  the  hot-wire  anemometer  system  accomplished  before 
and  after  data  collection  showed  that  the  changes  of  the  tunnel  mean 
velocity  as  a function  of  time  were  also  negligible. 


3.  UNSTEADY  FLOW  FIELD  HARMONIC  OPTIMIZATION 


Examination  of  the  unsteady  boundary  layer  equations  shows 
the  strong  influence  of  higher  harmonics  in  the  free  stream 
velocity  on  the  unsteady  flow  effects.  A detailed  presentation 
and  discussion  of  these  equations  is  given  in  section  4.2.  The 
goal  of  the  unsteady  flow  optimization  was  to  provide  a longitu- 
dinally oscillating  velocity  at  a desired  amplitude  ratio,  reduced 
frequency,  and  mean  velocity,  whose  power  spectrum  was  as  free  as 
possible  of  harmonic  content.  Of  the  four  controllable  tunnel 
parameters,  a combination  of  vane  size,  tunnel  steady  velocity,  and 
vane  rotation  frequency  (or  rotation  period)  were  selected  to  give 
the  desired  reduced  frequency  and  velocity  oscillation  amplitude 
i„tio.  The  upper  and  lower  bypass  openings,  each  at  a setting  of 
L/2,  were  then  adjusted  for  maximum  second  harmonic  suppression. 

The  variable  bypass  was  intended  to  provide  a more  uniform  flow  to 
the  fan  over  the  period  of  vane  rotation.  The  relationship  of  the 

bypass  resistance  to  the  tunnel  flow  resistance  is  described  in  the 
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flow  analysis  by  McMichael  and  summarized  in  Appendix  B.  This 
analysis  showed  that  the  smallest  temporal  change  in  the  pressure 
rise  across  the  fan  occurred  when  the  pressure  and  frictional  losses 
across  the  bypass  ports  and  vanes  were  equal  to  the  sum  of  the 
pressure  drop  across  the  main  tunnel  vanes  and  the  frictional  losses 
in  the  test  section.  Careful  adjustment  of  the  bypass  port  openings 
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permitted  a fine  tuning  of  the  bypass  frictional  losses  and 
permitted  an  almost  steady  flow  to  reach  the  fan. 

The  actual  unsteady  flow  characteristics  generated  by  the 
smallest  set  of  rotating  vanes  (set  A)  can  be  compared  most  easily 
to  the  analysis  because,  for  this  case,  the  size  of  the  tunnel  and 
bypass  vanes  were  the  same.  The  analysis  showed  that  this  arrange- 
ment provided  a great  simplification  and,  ideally,  generated  only 
even  harmonics  in  the  oscillating  velocity  field.  The  absence  of 
add  harmonics  was  difficult  to  confirm  experimentally  since,  for 
the  small  vanes,  usually  only  the  second  harmonic  was  large  enough 
to  be  detected  by  the  Fourier  Analyzer  on  a logrithmic-magnitude 
power  spectrum  (maximum  scale  of  80  db). 

A typical  power  spectrum  of  the  linearized  velocity  signal 
for  the  small  vanes  (set  A)  is  shown  in  Figure  21  for  k = 0.873  and 
L/2  = 1 inch  (2.54  cm).  The  velocity  data  had  been  ensemble- 
averaged  over  50  data  block  entries.  The  measurements  were 
conducted  in  an  empty  test  section,  without  the  wing  installed. 

The  second  harmonic  power  was  observed  to  be  approximately  48  db 
lower  than  the  power  of  the  fundamental.  The  db  differences  between 
the  voltage  power  levels  of  the  first  and  second  harmonics, 

20  log-|Q  (Vj/V2),  are  plotted  in  Figure  22  for  a variation  of  L/2 
at  this  flow  condition.  The  optimum  or  maximum  power  difference  is 
observed  to  be  approximately  48  db  and  the  smallest  difference  is 
32  db.  The  power  differences  between  the  fundamental  and  third 
harmonic  and  between  the  fundamental  and  fourth  harmonic  were  also 
noted,  if  measurable.  In  no  case  were  they  measured  to  be  less  than 
60  db  and  70  db  respectively  for  the  small  vanes. 


DIFFERENCE  IN  FIRST  AND  SECOND  HARMONIC  POWER  LEVELS 
Small  Vanes  (set  A),  k = 0.0873  (2  Hz) 
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The  power  spectrum  results  for  the  larger  vane  sets  were 
acceptable,  but  not  as  free  from  harmonics  as  those  for  the  smallest 
vanes.  The  power  differences  between  the  fundamental  and  second 
harmonic  as  functions  of  L/2  are  presented  in  Figure  23  for  the 
medium  vanes  (set  B)  and  a mean  velocity  of  31  feet  per  second 
(9.45  meters  per  second)  and  in  Figure  24  for  the  large  vanes  (set 
C)  and  a mean  velocity  of  12  feet  per  second  (3.66  meters  per 
second).  These  data  were  collected  with  the  airfoil  installed  in 
the  test  section.  At  these  flow  conditions,  the  medium  and  large 
vane  sets  were  observed  to  produce  maximum  harmonic  differences  of 
15  db  and  11  db,  respectively.  The  maxima  can  be  clearly  distin- 
guished, but  they  are  not  as  large  in  magnitude  as  those  measured 
with  the  small  vane  set  (set  A)  installed.  It  must  be  remembered 
that,  for  mechanical  simplicity,  only  the  main-flow  vanes  were 
replaced  when  changing  vane  sets.  The  bypass  vanes  were  not 
replaced.  This  arrangement  greatly  complicated  the  flow  analysis, 
and  the  relationships  developed  in  Appendix  B do  not  clearly  show 
that  a maximum  exists  for  the  second  harmonic  suppression.  The 
experimental  results  show  that,  even  with  different  sizes  of  main 
and  bypass  vanes,  the  maxima  can  be  achieved  with  the  proper  setting 
of  L/2. 

In  order  to  obtain  the  parameters  required  for  the  flow 
reversal  experiments,  the  optimum  setting  of  L/2  could  not  always 
be  used.  The  objectives  of  the  flow  reversal  experiments  required 
two  separate  parametric  variations  of  the  flow.  The  first  varied 
the  reduced  frequency  and  held  the  mean  velocity  and  oscillation 
amplitude  nearly  constant.  The  second  varied  the  oscillation 

L, j 


Figure  24 

FIRST  AND  SECOND  HARMONIC 
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amplitude,  with  the  reduced  frequency  and  mean  velocity  almost 
constant.  It  was  observed  that  the  magnitudes  of  the  mean  velocity 
and  amplitude  ratio  varied  considerably  with  L/2  for  a given 
frequency.  The  flow  parameters  observed  while  measuring  the  flow 
reversal  data  described  in  Chapter  4 are  shown  in  Table  II.  Even 
though  off-optimum  settings  of  L/2  were  used,  the  second  harmonic 
rejection  was  at  least  10  db  for  all  cases  except  the  lowest 
frequency  (k  = 0.119). 

TABLE  II 

FLOW  PARAMETERS  - FLOW  REVERSAL  MEASUREMENTS 


Vane  L/2  U0  4 

Set  in  ft/sec  (sec  I ft/sec  I % 


R 

xlO5 

Vx-V2 
db 

2.421 

47 

1.717 

49 

1.383 

18 

1.459 

16 

1.470 

10 

0.570 

13 

0.542 

11 

0.5  60 

45 


B 7.0  60 


3.0  40 

60 


As  an  aside,  it  was  found  that  harmonic  optimization  could 
also  be  confirmed  using  the  method  of  Charnay  and  Mathieu  ® who 
measured  the  skewness  and  kurtosis  of  the  longitudinal  turbulent 
velocity  component,  u '.  They  found  the  conditions  where  S 
approached  zero  and  K approached  l.j  represented  an  optimum. 
Figures  25  and  26  show  the  values  of  K and  S versus  L/2  with  the 
small  vanes  (set  A)  for  the  same  flow  conditions  as  considered 
above.  It  can  be  seen  that  the  optimum  values  of  L/2  are  the  same 
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as  predicted  by  the  spectral  analysis.  Also,  the  points  of  maximum 
S and  K correspond  to  the  points  of  minimum  harmonic  suppression. 
The  zero  value  of  S indicates  that  the  distribution  of  the  velocity 
harmonics  is  Gaussian  and  that  the  disturbance  caused  by  the  vanes 
only  affects  the  pressure  field.  The  low  value  of  K shows  that  the 
Uj*  distribution  is  broad  for  the  optimized  tunnel  condition. 


4.  RESULTS  AND  DISCUSSION  OF  FLOW  REVERSAL  MEASUREMENTS 

4.1  Steady  Flow  Boundary  Layer  Reversal 

Boundary  layer  flow  reversal  was  measured  for  a variety  of 

steady  flow  conditions  prior  to  examining  the  unsteady  case.  The 

steady  data  were  considered  to  be  the  zero-frequency  limit  of  the 

unsteady  flow  measurements.  These  regions  of  reversed  flow,  which 

occurred  near  the  airfoil  leading  edge,  represented  the  lower 

portion  of  the  so-called  laminar  separation  bubble.  A model  of  this 

bubble  is  presented  in  Figure  27  taken  from  reference  31.  This 

region  of  separated  flow  is  caused  by  the  large  curvature  of  the 

airfoil  surface  near  the  leading  edge.  After  the  laminar  boundary 

layer  separates,  it  transitions  to  turbulent  flow  and  then  reattaches 

to  the  airfoil  surface.  The  dividing  streamline  shown  in  Figure  2/ 

defines  the  boundaries  of  the  recirculating  region  contained  in  the 

bubble.  Flow  reversal  measurements  obtained  in  this  experiment 

determined  the  position  of  the  reversal  boundary  shown  in  the 

20 

figure.  Previous  data  for  the  NACA  0012  airfoil  have  shown  a 
separation  bubble  from  the  leading  edge  to  approximately  twenty 
percent  chord  at  an  angle  of  attack  of  11  degrees  and  a Reynolds 
number  of  250,000. 

A typical  steady  flow  reversal  observed  in  this  experiment  is 


shown  in  Figure  28.  The  vertical  axis  represents  the  non-dimensional 
height  of  the  sensor  above  the  airfoil  surface  and  the  horizontal 


Figure  28 

FLOW  REVERSAL  REGION 

Steady  Flow,  U = 31  ft/sec,  a = 10  degrees 
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axis  represents  the  non-dimensional  distance  of  the  sensor  from  the 
airfoil  leading  edge.  The  region  below  the  line  shown  on  the  graph 
represents  a region  of  reversed  flow.  Figure  29  shows  the  steady 
flow  reversal  for  other  angles  of  attack  at  the  same  free  stream 
velocity  of  31  feet  per  second  (9.45  meters  per  second).  These  data 
were  used  to  define  the  steady-flow  stall  angle  of  attack  (10.5 
degrees)  as  the  angle  where  the  flow  on  the  entire  airfoil  upper 
surface  becomes  reversed.  Figure  30  shows  similar  results  for  a 
lower  free  stream  velocity  of  12  feet  per  second  (3.66  meters  per 
second).  The  steady-flow  stall  angle  of  attack  was  measured  to  be 
unchanged  at  this  lower  velocity.  This  would  be  expected  since  the 
airfoil  stall  is  a function  only  of  angle  of  attack  and  is  insensi- 
tive to  the  small  change  in  Reynolds  number  between  these  two 
velocities. 

The  steady-flow  reversal  measurements  provided  a comparison 
for  the  subsequent  unsteady  flow  measurements.  The  unsteady  flow 
reversal,  described  below,  indicates  that  the  reversed  flow  region 
starts  near  the  leading  edge,  in  a manner  similar  to  the  steady  flow 
case,  and  then  spreads  downstream.  This  extended  flow  reversal 
feeds  mass  upstream  towards  the  leading  edge. 

4.2  Unsteady  Boundary  Layer  Equations 

The  effects  of  the  velocity  oscillation  frequency  and 
amplitude  on  the  reversal  regions  on  the  airfoil  may  be  more  clearly 
understood  by  considering  the  unsteady  boundary  layer  equations. 

33 

These  are  developed  below  following  the  discussion  of  McCroskey, 
who  considers  the  unsteady  terms  to  be  perturbations  to  a 
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quasi-steady  flow.  A quasi-steady  boundary  layer  response  is  one 
which  is  always  in  phase  with  the  oscillating  free  stream  velocity. 
The  equations  developed  below  were  later  non-dimensional ized  to 
observe  the  relative  influence  of  the  free  stream  velocity 
oscillation  frequency  and  amplitude.  The  unsteady  boundary  layer 
momentum  equation  is: 
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The  equation  for  the  undisturbed  free  stream  rnay  be  written 
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= 0. 


Substituting  this  into  equation  (1)  gives: 
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(2) 


McCroskey  interprets  the  first  term  on  the  right  hand  side  as  the 

contribution  of  the  quasi-steady  pressure  gradient.  The  second 

term  represents  the  unsteady  temporal  pressure  gradient  and  the 

third  term  is  interpreted  as  the  unsteady  temporal  viscous  effect. 

McCroskey  also  notes  that  the  second  and  third  terms  may  be 

considered  as  auxiliary  spatial  pressure  gradients  caused  by  the 

flow  unsteadiness.  This  interpretation  of  the  unsteady  terms  must 

be  modified  by  consideration  of  the  boundary  layer  geometry.  The 

influence  of  the  so-called  unsteady  temporal  viscous  effect  is 

important  only  at  the  airfoil  surface  and  in  the  thin  viscous  shear 

2 

layer  described  above  by  Lighthill.  The  boundary  layer  region 
above  this  shear  layer  will  be  more  strongly  affected  by  the 
temporal  and  spatial  gradients.  Also,  the  results  described  above 


j 
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obtained  by  Simpson15  have  indicated  that  the  effects  of  turbulent 
mixing  must  be  considered  after  boundary  layer  separation  has 
occurred.  The  boundary  layer  equation  resulting  from  including  the 
turbulent  velocities  would  include  additional  unsteady  terms  due  to 
the  Reynolds  stresses. 

For  the  moment,  the  external  oscillating  flow  velocity  is 
assumed  to  be 


U„(x,t)  = Uj  (x)  [1  + N sin  (out)]  (3) 

where  Uj  is  the  steady  mean  velocity  and  N is  the  amplitude  ratio, 
AU/Uj.  This  oscillating  velocity  may  be  used  to  determine  the 
effect  of  the  free  stream  oscillations  on  equation  (2).  After 
substitution,  the  second  term  on  the  right  hand  side  of  equation  (2) 
becomes : 

= Uj  (x)  N w cos  (wt)  = Uj  (x)  N co  sin  (wt  + ~/2) 

= Aj (x)  sin  (wt  + n/2) . (4) 

The  unsteady  temporal  pressure  gradient  is  observed  to  lead  the 
phase  of  the  external  velocity  by  a quarter-period.  The  first  term 
on  the  right  hand  side  of  equation  (2)  becomes: 
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dx 


sin7  (wt) 


= Bj(x)  + Cj(x)  sin  (wt)  + Dj(x)  sin7  (wt) . (5) 


Substitution  of  equations  (4)  and  (5)  into  equation  (2),  using  the 
following  non-dimensionalization,  yields: 
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3u‘  , a2u 1 _ 9...  9,  9u 1 . l du, 

U 3x‘  V ^v/'2  2kN  cos  (t  ) - 2k  , + |,  Axt% 


3y  3y 


at ' u,  dx' 


. 2N  dUi  . 4 ,,  , N2  dUj  . 2 l4.,\ 

+ u;  d^51"  (t  } + u;dx4sin  (t  } 


= A2  sin  (t 1 + tt/2)  - 2k  |p-  + B2(x) 


+ C (x)  sin  (t1)  + D (x)  sin2  (t‘) 


(6) 


where: 


The  first  two  terms  on  the  right  hand  side  of  equation  (6)  represent, 
in  a non-dimensional  form,  the  unsteady  temporal  gradient,  and  the 
last  three  terms  represent  the  unsteady  spatial  gradient.  The 
airfoil  geometry  provides  a boundary  condition  which  determines 
whether  the  unsteady  spatial  pressure  gradient  is  favorable 
(negative)  or  adverse  (positive).  The  airfoil  angle  of  attack 
provides  a boundary  condition  which  determines  the  relative  magni- 
tude of  the  favorable  or  adverse  condition.  The  favorable  or 
adverse  nature  of  the  spatial  unsteady  pressure  gradient  depends 
principally  on  the  sign  of  the  coefficient  C2(x),  for  small  values 
of  the  amplitude  ratio,  N.  For  large  oscillation  amplitude  ratios 
near  one,  however,  the  coefficient  D?(x)  becomes  equally  important. 
Large  oscillation  amplitude  ratios  can  have  a large  influence  on  the 
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phase  lead  or  lag  of  the  unsteady  boundary  layer  response  through 
the  trigonometric  relationship: 

D2(x)  sin2V  = [i  + cos  2t'J. 

The  three  terms  of  the  unsteady  spatial  gradient  are  observed 
to  be  functions  of  the  mean  velocity,  Uj,  and  velocity  oscillation 
amplitude  ratio,  N,  but  not  of  the  reduced  frequency,  k.  The 
unsteady  temporal  pressure  gradient,  which  is  the  first  term  on  the 
right  hand  side  of  equation  (6),  depends  on  both  the  reduced 
frequency  and  the  oscillation  amplitude  ratio.  Because  of  differen- 
tiation, the  contribution  of  this  term  leads  the  phase  of  the  free 
stream  velocity  by  tt/2  radians.  The  influence  of  this  term  is 
opposed  by  the  effect  of  the  unsteady  viscous  term,  the  second  term 
on  the  right  hand  side  of  equation  (6).  The  viscous  term  is 
observed  to  be  a function  of  the  reduced  frequency,  and  is  indepen- 
dent of  the  oscillation  amplitude  ratio.  This  term  will  be  of  very 
small  influence  away  from  the  airfoil  surface. 

The  effect  of  the  unsteady  spatial  and  temporal  gradients  on 

32 

the  shear  at  the  wall  are  summarized  in  Figure  31,  from  McCroskey. 
The  graphs  on  the  left  hand  side  of  the  figure  show  the  influence 
of  a favorable  unsteady  spatial  pressure  gradient  (C?(x)  >0),  while 
the  right  hand  sides  show  the  effect  of  an  adverse  gradient 
(C2(x)  < 0).  This  unsteady  contribution  is  similar  to  the  effect 
of  an  adverse  spatial  pressure  gradient  in  steady  flow,  which  causes 
steady  boundary  layer  separation.  The  unsteady  adverse  pressure 
gradient,  however,  is  shown  to  be  of  opposite  phase  as  well  as 
opposite  sign.  The  graphs  of  wall  shear  at  the  bottom  of  Figure  31 
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indicate  how  the  unsteady  gradients  cause  the  shear  to  depart  from 
the  quasi -steady  response.  This  presentation  shows  that  the  shear 
leads  the  phase  of  the  quasi-steady  response  for  a favorable, 
unsteady  spatial  pressure  gradient,  and  lags  the  quasi-steady 
response  for  an  adverse  spatial  gradient.  It  must  be  remembered 
that  this  shear  was  produced  by  a small  oscillation  amplitude  ratio. 

The  actual  phase  lead  or  lag  will  ultimately  depend  on  whether  the 
spatial  or  temporal  gradients  dominate  the  boundary  layer  response. 

Both  phase  leads  and  lags  from  the  quasi-steady  response  were  noted 
in  the  experimental  data  presented  below.  These  data  use  flow 
reversal  in  the  velocity  field  measured  above  the  airfoil  surface, 
rather  than  wall  shear,  to  indicate  the  onset  of  unsteady  flow 
effects. 

The  manner  in  which  the  term  "quasi -steady"  is  used  in  this 

experiment  should  be  more  clearly  defined,  since  the  term  has 

several  different  meanings  in  the  literature.  Some  authors^ 

consider  a quasi-steady  boundary  layer  response  to  be  one  which  is 

1 7 

in  phase  with  the  time-varying  external  velocity.  Others  define 
a quasi-steady  response  as  one  which  is  identical  in  magnitude  to  an 
equivalent  series  of  steady-flow  responses  measured  for  the  fixed 
values  of  the  free  stream  velocity  occurring  at  each  phase  of  the 
unsteady  velocity  oscillations.  In  this  experiment,  for  flow 
reversal  measured  at  the  airfoil  surface,  the  flow  reversal  is 
defined  to  be  quasi-steady  if  it  initially  appears  when  the 
instantaneous  free  stream  velocity  magnitude  drops  below  the  mean 
velocity  magnitude,  and  disappears  when  the  velocity  magnitude 
increases  above  the  mean.  If  the  onset  of  flow  reversal  is  noted 

k - 
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before  or  after  the  phase  where  the  varying  free  stream  velocity 
passes  through  the  mean  value  (t‘  = 0,tt)  , it  is  said  to  lead  or  lag 
the  phase  of  the  free  stream  velocity.  A quasi-steady  flow  reversal 
response  is  more  difficult  to  describe  for  regions  above  the  airfoil 
surface.  The  reversal  region  is  said  to  be  quasi-steady  if  the 
region  is  symmetrical  with  respect  to  the  free  stream  velocity.  The 
reversal  region  will  have  a maximum  vertical  height  when  the  magni- 
tude of  the  free  stream  velocity  is  a minimum.  If  the  reversal 
region  is  skewed  to  the  left  or  right,  it  is  said  to  lead  or  lag, 
respectively,  the  phase  of  the  free  stream. 

If  the  assumed  form  of  the  free  stream  velocity  given  in 
equation  (3)  is  not  correct,  that  is,  if  the  flow  is  not  sinusoidal 
but  varies  as: 


Ujx.t)  = Uj  (x) 


1 + N 


|an  cos  (nmt)  + bn  sin  (neat) 


(7) 


then  the  momentum  equation  must  be  altered  as  follows: 

oo 

u'  fxT  + v'  fyr  ■ = 2kN  2 [nbn  Sin  (nt‘  + tt/2)  - na()  sin  (nt')J 
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It  can  be  seen  in  equation  (8)  that  the  viscous  term,  the  second 
term  on  the  right  hand  side,  is  not  affected  by  the  higher 
harmonics.  The  higher  harmonics  manifest  themselves  as  pressure 
waves  and  not  as  viscous  shear.  The  higher  harmonic  terms  could, 
however,  drastically  alter  the  influence  of  both  the  spatial  and 
temporal  unsteady  terms.  The  complicated  interaction  of  the  high 
frequency  terms  in  equation  (8)  shows  the  desirability  of  having  the 
free  stream  velocity  oscillations  as  free  as  possible  of  high 
harmonic  content. 

A more  detailed  investigation  should  consider  the  oscillating 
velocity  at  the  outer  edge  of  the  airfoil  boundary  layer,  of  the 
form 

Uoo(x,t)  = Uj (x)  1 1 + M ( x ) sin  [wt  + G(x)]( 

where  both  the  amplitude  ratio,  N,  and  a phase  angle,  0,  are 
functions  of  distance  along  the  airfoil  chord.  After  substitution, 
several  additional  terms  are  observed  in  the  non-dimensional 
boundary  layer  equation: 


cos  (t'  +0)  + N2  f-,  sin  (t*  + 0)  cos  (t1  +0). 
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Using  the  Fourier  series  presentation: 


Ujx.t)  = Ut(x) 


1 + N(x)  X]  [an(x)  cos  (n^t)  + bp(x)  sin  (ncat)l 
n-1  *-  J 


the  non-dimensional  boundary  layer  equation  becomes 
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The  additional  terms  described  above  were  considered  to  be 
very  small  since  the  wavelength  of  the  free  stream  velocity 
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oscillation  was  longer  than  the  length  of  the  wind  tunnel  test 
section.  Preliminary  correlation  measurements  of  velocity  near  the 
airfoil  leading  and  trailing  edges  showed  almost  exact  correlation. 

4.3 Unsteady  Flow  Reversal  for  Angles  of  At tack  Below  th e Steady 

Flow  Stall  Angle 

The  locus  of  points  which  indicated  the  onset  of  reversed 
flow  was  used  to  define  the  flow  reversal  region.  Figure  32 
presents  unsteady  flow  reversal  measurements  for  an  airfoil  angle 
of  attack  of  9 degrees  (below  the  measured  steady  flow  stall  angle 
of  10.5  degrees)  and  a mean  velocity  of  31  feet  per  second  (9.45 
meters  per  second).  The  spatial  extent  of  the  flow  reversal  region, 
expressed  by  nondimensional  chord  and  height  above  the  airfoil 
surface,  is  indicated  for  several  phase  points  in  the  oscillation 
period.  The  regions  below  the  curves  represent  reversed  flow.  The 
limited  upstream  extent  of  the  reversal  measurements  was  due  to  the 
geometry  of  the  split-film  probe.  Although  data  was  collected  for 
256  phase  points  in  the  half-period  shown,  the  limitations  of  the 
graphical  data  presentation  dictated  that  only  a small  number  of 
phase  points  be  plotted.  Two-dimensional  presentations  of  these 
data,  shown  below,  include  data  for  more  phase  points.  Figure  33 
presents  similar  data  collected  at  a higher  reduced  frequency,  but 
with  approximately  the  same  oscillation  amplitude  ratio.  Flow 
reversal  measured  at  this  same  angle  of  attack,  but  at  a much  larger 
oscillation  amplitude  ratio,  are  shown  in  Figure  34. 

For  this  angle  of  attack,  the  reversed  flow  region  was 
observed  to  be  similar  in  size  and  shape  to  the  flow  reversal  which 
formed  the  separation  bubble  noted  in  steady  flow.  The  flow 
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reversal  was  observed  only  during  the  second  half  of  the  velocity 
oscillation  period  (-n  £ t'  £ 2tt)  when  the  magnitude  of  the  free 
stream  velocity  was  smaller  than  the  mean  velocity.  At  other 
oscillation  phase  points  the  reversed  flow  disappeared,  indicating 
that  the  laminar  flow  near  the  airfoil  leading  edge  never  separated. 

4.4  Unsteady  Flow  Reversal  for  Angles  of  Attack  Above  the  Steady 
Flow  Stall  Angle 

Reversal  data  measured  at  the  higher  angles  of  attack  indi- 
cated that  the  flow  reversal  region  extended  over  the  entire  airfoil 
upper  surface  and,  in  many  cases,  thickened  near  the  trailing  edge. 
The  adverse  unsteady  pressure  gradient  appeared  to  exert  a stronger 
influence  on  the  onset  of  flow  reversal  than  noted  at  lower  angles 
of  attack.  Measurements  were  made  at  an  angle  of  attack  of  11 
degrees,  just  above  the  steady  flow  stall  angle,  and  also  at  an 
angle  of  attack  of  12  degrees.  Reversal  data  shown  in  Figure  35 
were  mtasured  at  an  intermediate  reduced  frequency  (vane  rotation 
speed  of  2.0  Hz)  and  a low  oscillation  amplitude  ratio.  Similar 
data  for  two  higher  frequencies  are  presented  in  Figures  36  and  37. 
Figure  38  shows  measurements  made  at  a much  lower  reduced  frequency 
of  k = 0.119,  produced  by  using  r.  lower  vane  rotation  speed  of  1.43 
Hz.  An  increase  in  the  angle  of  attack  to  12  degrees  produced  the 
flow  reversals  presented  in  Figures  39  and  40,  measured  at  inter- 
mediate and  high  reduced  frequencies,  respectively.  The  influence 
of  a large  oscillation  amplitude  ratio  was  measured  only  at  an 
angle  of  attack  of  11  degrees  and  is  shown  in  Figure  41. 

It  was  observed  that  the  flow  reversal  region  initially 
appeared  near  the  airfoil  leading  edge.  This  result  is  in  contrast 


Figure  35 
low  rever: 

legrees,  k 


ire  36 
REVER! 


ire  38 

rever: 


Figure  39 
Y FLOW  REVER: 
2 degrees,  k 


Figure  41 

UNSTEADY  FLOW  REVERSAL  REGION 
11  degrees,  k = 0.168,  N = 0 


85 


to  trailing  edge  flow  reversal  noted  by  Carr,  McAlister  and 
McCroskey3  for  the  oscillating  NACA  0012  airfoil.  At  first 
consideration,  the  leading  edge  flow  reversal  would  appear  to  be 
closer  in  nature  to  the  results  noted  in  Figure  2 for  the  NACA  0012 
airfoil  with  the  sharp  leading  edge.  However,  this  leading  edge 
reversal  phenomena  was  identified  as  a "bursting"  of  the  leading 
edge  laminar  bubble,  and  was  not  considered  the  classical  dynamic 
stall  flow  model.  It  is  felt  that,  for  the  stationary  airfoil  in 
the  oscillating  flow,  the  initial  onset  and  downstream  growth  of 
the  flow  reversal  region  is  due  to  the  energy  being  transferred  to 
the  laminar  reversal  region  by  the  turbulent  flow.  The  reversal 
region  is  not  due  to  a downstream  extension  of  the  leading  edge 
laminar  bubble.  It  is  also  felt  that  the  trailing  edge  reversal 
was  not  observed  on  the  stationary  airfoil  because  of  the  absence 
of  pitch  rate  effects.  The  increased  unsteady  spatial  pressure 
gradients  near  the  trailing  edge  caused  by  the  inertia  of  the 
airflow  close  to  the  moving  airfoil  were  not  present. 

The  shape  of  the  reversal  regions  described  in  the  figures 
above  bears  a strong  resemblance  to  the  reversal  and  separation 
modH  of  Sears  and  Telionis,^  shown  in  Figure  3.  The  thickening 
of  the  flow  reversal  near  the  trailing  edge  in  Figure  37  may  be 
indicative  of  unsteady  boundary  layer  separation.  The  larger 
reversal  region  for  the  larger  angle  of  attack  of  12  degrees  would 
seem  to  indicate  that  unsteady  separation  has  occurred  further 
upstream.  The  model  described  in  Figure  37  considered  only  one 
phase  point  of  the  unsteady  flow  period  and  did  not  include  the 
formation  of  any  leading  edge  vortex  disturbances.  Higher 


frequencies  were  observed,  in  Figure  39,  to  eliminate  the  reversal 
region  thickening  near  the  trailing  edge.  According  to  the  above 
model,  this  would  represent  the  absence  of  unsteady  boundary  layer 
separation  throughout  the  oscillation  period.  This  is  an  important 
result  when  applied  to  the  problem  of  dynamic  stall.  At  higher 
oscillation  frequencies,  the  flow  reversal  could  be  responsible  for 
feeding  mass  to  the  leading  edge  vortex-like  disturbance  without 
separation  occurring  near  the  trailing  edge. 


5.  FREQUENCY  EFFECTS  ON  FLOW  REVERSAL 
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5.1  Frequency  Effects  for  Angles  of  Attack  Below  the  Steady  Flow 
Stall  Angle 

The  changes  of  the  shape  of  the  flow  reversal  region 
produced  by  varying  the  reduced  frequency  of  oscillation  are  indi- 
cated in  Figure  42  for  an  airfoil  angle  of  attack  of  9 degrees. 

Close  to  the  airfoil  surface,  there  appeared  to  be  very  little 
influence  caused  by  the  flow  unsteadiness,  since  the  chordwise 
extent  of  the  flow  reversal  was  almost  the  same  for  all  frequencies, 
including  steady  flow.  As  the  reduced  frequency  was  increased, 
however,  the  reversal  region  increased  in  height  above  the  airfoil 
surface. 

These  height  increases  are  believed  to  be  caused  by  changes 
in  the  unsteady  centrifugal  forces  along  the  streamlines,  rather 
than  the  unsteady  terms  described  in  equation  (6).  The  unsteady 
effects  may  be  interpreted  using  boundary  layer  equations  written 
in  generalized  orthogonal  coordinates  aligned  with  the  airfoil 
surface.  Consider  first  the  steady  momentum  equation  for  flow 
normal  to  the  airfoil  surface  (along  the  y axis): 
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(9) 


This  equation  indicates  the  pressure  gradient,  on  the  right  hand 
side,  necessary  to  balance  the  centrifugal  force,  on  the  left  hand 
side.  These  forces  are  due  to  the  local  streamline  curvature,  k, 
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above  the  airfoil  surface.  When  equation  (9)  is  satisfied,  the 
steady  flow  follows  the  curvature  of  the  airfoil  and  is  not 
reversed.  The  steady  flow  data  in  Figure  42  show  that  the  flow 
was  reversed  ahead  of  the  chord  position  x = 0.14.  For  this  small 
region,  the  centrifugal  forces  were  greater  than  the  pressure 
gradient  and  equation  (9)  was  not  satisfied.  It  appears  that  a 
separation  bubble  was  formed,  as  indicated  by  the  presence  of  the 
reversed  flow  near  the  leading  edge  for  steady  flow  conditions. 

The  chordwise  extent  of  the  bubble  is  limited  both  by  a decrease 
in  the  local  airfoil  curvature  and  by  a transition  to  turbulent 
flow.  Experimental  data  showed  the  boundary  layer  to  be  turbulent 
downstream  of  the  reattachment  point. 

For  unsteady  flow,  another  term  appears  in  the  normal-flow 


momentum  equation: 
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(10) 


The  local  acceleration  term  augments  the  centrifugal  force  that 
must  be  balanced  by  the  pressure  gradient.  Higher  oscillation 
frequencies,  acting  through  this  unsteady  term,  will  produce  larger 
imbalances  in  equation  (10)  and  larger  flow  reversal  regions.  The 
high  frequency  data  in  Figure  30  confirm  this  trend.  These  data 
also  show  the  strong  influence  of  the  wall  curvature.  For  the 
lower  curvatures  occurring  aft  of  x = 0.17,  no  flow  reversals 
occurred  at  any  reduced  frequency. 

Figure  43  shows  the  change  in  phase  of  the  reversed  flow 
region  over  one  period  of  the  velocity  oscillation.  Reversal 
height  data  for  the  two  frequencies  were  plotted  for  a fixed 


chord-wise  position  of  10  percent.  The  phase  relationships  of  the 
data  at  both  frequencies  appeared  to  be  close  to  quasi -steady , 
that  is,  in  phase  with  the  oscillating  free  stream.  Both  curves 
were  skewed  slightly  to  the  left,  due  to  the  effect  of  the  favorable 
unsteady  spatial  pressure  gradients  near  the  airfoil  leading  edge. 
The  higher  reduced  frequency  increased  the  phase  lead,  but  not 
through  the  unsteady  temporal  gradients.  This,  again,  is  most 
probably  an  effect  of  the  unsteady  contribution  to  the  centrifugal 
force  along  the  streamline. 

5,2  Frequency  Effects  for  Angles  of  Attack  Above  the  Steady  Stall 
Angl  e 

The  influence  of  changes  in  the  reduced  frequency  on  the 
shape  of  the  flow  reversal  region  is  observed  in  Figure  44.  Data 
are  presented  for  an  angle  of  attack  of  11  degrees  and  for  several 
phase  points  in  the  oscillation  period.  For  the  low  frequency  case 
shown,  the  flow  reversal  region  extended  over  the  entire  airfoil 
surface  and  the  reversal  region  increased  in  height  near  the 
trailing  edge.  The  region  of  flow  reversal  ahead  of  the  mid-chord 
position  grew  due  to  the  reversed  flow  from  the  trailing  edge.  The 
spatial  extent  of  the  reversal  region  to  the  airfoil  trailing  edge 
indicates  that  the  adverse  unsteady  spatial  gradients  had  a 
dominant  effect  on  the  flow. 

The  data  presented  in  Figure  44  for  the  higher  frequency 
show  a reduction  in  size  of  the  flow  reversal  region.  The  flow 
reversal  near  the  trailing  edge  was  entirely  eliminated.  The 
unsteady  temporal  gradients  from  equation  (6)  were  observed  to  be 
directly  proportional  to  reduced  frequency  and  thus  made  a larger 


contribution  to  the  unsteady  flow  response  at  higher  frequencies, 
opposing  the  spatial  gradients. 
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The  shape  of  reversal  data  collected  at  an  angle  of  attack 
of  12  degrees  for  two  reduced  frequencies  are  compared  in  Figure  45. 
The  influence  of  the  high  reduced  frequency  is  not  as  pronounced  at 
this  larger  angle  of  attack  since  the  unsteady  boundary  layer  is 
separated  for  part  of  the  oscillation  cycle.  While  the  maximum 
heights  of  the  reversal  regions  for  the  two  frequencies  are  closer 
in  magnitude,  the  higher  frequency  effects  are  again  observed  to 
limit  the  chordwise  extent  of  the  reversal  region.  The  important 
influence  of  the  angle  of  attack  increase  is  observed  with  the  phase 
of  the  reversal  region  response  and  is  shown  below. 

The  influence  on  the  phase  due  to  changes  in  reduced 
frequency  was  observed  to  be  more  significant  above  the  steady  stall 
angle.  The  reversal  data  are  compared  in  Figure  46  for  an  angle  of 
attack  of  11  degrees  and  several  chordwise  positions.  The  phase  of 
the  lower  frequency  reversed  flow  region  is  observed  to  lag  the 
quasi-steady  response  near  the  leading  edge,  but  is  close  to  the 
quasi-steady  phase  near  the  trailing  edge.  Near  the  leading  edge, 
the  influence  of  the  reduced  frequency,  acting  through  the  unsteady 
temporal  gradients,  appears  to  be  dominant  and  causes  the  phase  lag. 
The  effects  of  the  spatial  gradients  are  insignificant  since,  at 
this  forward  chord  position  and  low  amplitude  ratio,  the  spatial 
gradient  terms  are  small,  either  slightly  negative  (favorable)  or 
slightly  positive  (adverse).  Near  the  airfoil  trailing  edge,  the 
spatial  gradient  terms  are  larger  and  oppose  the  frequency-induced 
phase  lag.  The  higher-frequency  data  in  Figure  46  are  shown  to  have 
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*'  Figure  46 

FREQUENCY  EFFECT  ON  PHASE  OF  FLOW  REVERSAL  RESPONSE,  a = 11  degrees 


a phase  lag  at  all  chord  positions  except  at  the  trailing  edge. 

For  this  frequency,  the  temporal  gradients  are  the  dominant 
influence  on  the  reversal  region  response. 

An  increase  in  the  angle  of  attack  was  observed  to  produce 
important  changes  in  the  phase  of  the  reversed  flow  region.  Figure 
47  compares  the  phase  of  reversal  data  at  an  angle  of  attack  of  12 
degrees  for  two  reduced  frequencies.  For  both  frequencies,  the 
response  is  seen  to  lead,  rather  than  lag,  the  quasi-steady  phase. 
This  result  agrees  with  the  analysis  for  chord  positions  near  the 
leading  edge  where  the  unsteady  spatial  gradient  had  a favorable 
influence.  The  higher  reduced  frequency  data  had  a smaller  phase 
lead  due  to  the  temporal  gradient  efforts.  However,  the  flow 
reversal  data  for  the  aft  portion  of  the  airfoil  show  a phase  lead, 
even  though  the  spatial  gradient  influence  near  the  airfoil  surface 
had  become  adverse.  This  result  was  not  predicted  by  the  unsteady 
boundary  layer  equations  based  on  attached  flow. 

At  this  larger  angle  of  attack,  the  boundary  layer  was 
turbulent  and  probably  separated,  and  was  not  governed  by  the 
generalized  boundary  layer  equations  presented  above.  The  phase 
response  for  the  portion  of  the  flow  reversal  region  av/ay  from  the 
surface  led  the  free  stream  due  to  the  rapid  vertical  transfer  of 
energy  at  this  height.  Closer  to  the  surface,  the  flow  reversal 
tended  to  lag  the  free  stream  because  of  slower  energy  transfer 
because  of  viscous  effects  near  the  surface. 

Experiments  on  dynamic  stall  using  the  oscillating  airfoil 
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by  Carr,  McAlister,  and  McCroskey  indicated  that  an  increase  in 
the  reduced  frequency  of  oscillation  delayed  the  onset  of  flow 
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reversal  and  the  subsequent  series  of  dynamic  stall  events.  The 
strong  influence  of  high  frequency  on  the  trailing  edge  flow 
reversal  appears  to  be  present  for  both  the  oscillating  airfoil  and 
oscillating  free  stream  experiments.  Proper  control  of  this 
unsteady  parameter  may  provide  part  of  the  solution  to  the 
constructive  use  of  dynamic  stall  for  maneuvering  flight. 

A curious  result  is  observed  at  an  angle  of  attack  of  11 
degrees  for  the  lower  reduced  frequency  of  k = 0.119.  Data 
presented  earlier  in  Figure  38  showed  that  flow  over  the  entire 
airfoil  surface  was  reversed  for  much  of  the  oscillation  period. 

This  would  appear  to  indicate  a stronger  influence  of  the  unsteady 
spatial  gradients.  The  phase  relationships  plotted  in  Figures  48 
and  49  indicated  that  the  flow  reversal  leads,  rather  than  lags,  the 
quasi -steady  response  for  this  lower  frequency.  This  phase  lead  is 
contrary  to  the  phase  lag  that  would  be  expected  for  a strong 
adverse  unsteady  spatial  pressure  gradient.  It  is  felt  that  this 
phase  lead  is  due  to  the  strong  presence  of  higher  harmonics  in  the 
free  stream  velocity  oscillations  at  this  lower  frequency.  The 
higher  harmonics  are  seen,  in  equation  (8),  to  influence  both  the 
temporal  and  spatial  unsteady  pressure  gradients.  The  phase  lead  is 
due  to  the  higher  harmonic  coefficients  contained  in  the  unsteady 
temporal  pressure  gradient.  This  contribution  of  this  term  is  seen 
to  lead  the  quasi-steady  response  through  the  phase  factor  of  m/2. 
The  strong  influence  of  the  higher  harmonic  coefficients  on  the 
unsteady  spatial  pressure  gradient  terms  in  equation  (8)  is  seen  by 
the  increased  extent  of  reversed  flow  towards  the  airfoil  trailing 
edge. 


Ngure  49 

FREQUENCY  EFFECT  ON  PHASE  OF  FLOW  REVERS, 


These  low  frequency  results  emphasize  the  desirability  of 
conducting  unsteady  experiments  of  this  type  in  oscillating  flows 
that  are  free  from  higher  harmonic  content.  These  harmonics  can 
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produce  data  which  indicate  phase  leads  where  lags  are  normally 
observed.  Also  the  relative  influence  of  the  spatial  and  temporal 
unsteady  effects  can  be  greatly  altered  by  the  presence  of  the 
higher  harmonics. 


A 


6.  OSCILLATION  AMPLITUDE  EFFECTS  ON  FLOW  REVERSAL 


6.1  Amplitude  Effects  for  Angles  of  Attack  Below  the  Steady  Flow 
Stall  Angle 

Flow  reversal  data  presented  below  show  the  effect  on  the 
flow  reversal  caused  by  a large  increase  in  the  amplitude  ratio,  N, 
of  the  oscillating  free  stream.  The  influence  of  this  parameter  is 
significant  due  to  the  coupling  of  the  temporal  and  spatial  unsteady 
gradients  through  N.  This  influence  is  shown  in  equation  (6),  which 
is  repeated  below: 


If  ♦ if  - If  - 2kN  <*'>  - 2k  If  + ^ Sr 

+ B7^s1"  (t')+D7^'sinJ  (t'> 


= fl2  sin  (t'  + »/2)  - 2k  |jr  + B2(x) 


+ C2(x)  sin  (f)  + D2(x)  sin9  (t',.  (6) 

It  can  be  noted  that  a combination  of  an  amplitude  ratio  near  unity 
and  a low  mean  free  stream  velocity,  U19  will  cause  the  unsteady 
pressure  gradient  terms  (the  last  three  terms  on  the  right  hand 
side)  to  dominate  the  boundary  layer  response. 

This  lower  speed  flow  situation  provided  the  data  shown  in 
Figure  50  at  an  angle  of  attack  of  9 degrees  for  a fixed  phase  point 
of  t'  = 3tt/4  and  a mean  velocity  of  12  feet  per  second  (3.66  meters 
per  second).  This  mean  velocity  was  the  maximum  attainable  with  the 


Figure  50 

AMPLITUDE  RATIO  EFFECT  ON  SHAPE  OF  REVERSAL  REGION 
a = 9 degrees,  t'  = 3tt/2,  Uj  = 12  ft/sec 
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large  set  of  rotating  vanes  (set  C).  The  mean  velocity  obtained 
with  the  medium  size  vanes  (set  B)  was  then  adjusted  to  12  feet  per 
second.  It  is  observed  that,  for  the  larger  amplitude  ratio, 
reversed  flow  extends  to  the  trailing  edge  and  appears  to  be 
separating  near  the  trailing  edge.  For  this  flow  condition,  the 
effect  of  the  unsteady  adverse  pressure  gradient  appears  to  dominate 
the  reversed  flow  over  the  entire  airfoil  upper  surface. 

The  effect  of  the  amplitude  ratio  variation  on  the  phase  of 
the  flow  reversal  is  shown  in  Figure  51  for  a chord  position  of 
x = 0.15.  The  phase  lead  caused  by  the  large  amplitude  ratio  is 
significant  and  is  probably  due  to  the  higher  order  unsteady  spatial 
gradient  in  equation  (6): 

0 (x)  sin  (t')  = [1  +cos  (2t')]. 

For  the  flow  reversal  caused  by  the  lower  amplitude  ratio,  the  phase 
lead  is  due  to  the  unsteady  temporal  gradient. 

6.2  Amplitude  Effects  for  Angles  of  Attack  Above  the  Steady  Flow 
Stall  Angle 

The  effect  of  an  increase  in  oscillation  amplitude  ratio  on 
the  flow  reversal  region  shape  at  an  angle  of  attack  of  11  degrees 
was  not  as  pronounced  as  for  the  lower  angle  of  attack.  Figure  52 
shows  that,  for  a phase  point  of  t'  = 3tt/2,  the  reversal  regions 
produced  by  the  two  different  amplitude  ratios  did  not  differ 
greatly.  The  reversed  flow  region  extends  to  the  trailing  edge  for 
both  cases,  and  is  actually  thicker  for  the  lower  amplitude  ratio. 
This  increased  thickness  is  probably  due  to  the  presence  of  a larger 
region  of  separated  flow  for  the  lower  amplitude  ratio  case. 


Figure  5Z 

AMPLITUDE  RATIO  EFFECT  ON  SHAPE  OF  REVERSAL  RE 
a = 11  degrees,  t'  = 3tt/2,  U x = 12  ft/sec 
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For  this  angle  of  attack,  the  increased  amplitude  variation 
had  a smaller  effect  on  the  phase  of  the  flow  reversal.  Data  shown 
in  Figure  53  for  a fixed  chord  position  of  \ = 0.15  indicate  that 
unsteady  spatial  pressure  gradients  caused  a phase  lag  of  similar 
proportions  for  both  amplitude  ratios.  The  high  angle  of  attack, 
coupled  with  the  low  free  stream  velocity,  has  probably  caused  a 
flow  separation  during  the  second  half  of  the  oscillation  period. 

The  amount  of  vorticity  transferred  to  the  reversed  separated  flow 
appears  to  be  not  strongly  affected  by  the  amplitude  ratio,  but  the 
vorticity  transfer  extends  further  down  into  the  boundary  layer  for 


the  oscillations  with  the  larger  amplitude  ratio. 

Experimental  data  for  a variation  of  oscillation  amplitude 

for  the  oscillating  airfoil  were  measured  by  Carr,  McAlister,  and 
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McCroskey.  It  was  observed  that  an  increase  in  the  magnitude  of 
the  angle  of  attack  variation  caused  a phase  lead  in  the  onset  of 
flow  re /ersal  and  the  vortex  shedding  of  dynamic  stall.  This  phase 
lead  was  probably  caused  by  inertial  effects  caused  by  increases  in 
the  pitch  rate  of  the  airfoil,  and  are  not  directly  related  to  the 
amplitude  variations  considered  for  the  oscillating  free  stream. 


Figure  53 

AMPLITUDE  RATIO  EFFECT  ON  PHASE  OF  REVERSAL  REGION 
a = 11  degrees,  x = 0.15 


7.  CONCLUSIONS 

The  onset  of  unsteady  boundary  layer  flow  reversal  was 
measured  above  the  surface  of  a stationary  NACA  0012  airfoil  in  a 
longitudinally  oscillating  flow.  The  response  of  the  flow  reversal 
region  to  the  influence  of  the  flow  unsteadiness  was  found  to  be 
similar  to  the  flow  reversal  observed  during  dynamic  stall. 

Previous  dynamic  stall  experiments  using  oscillating  airfoils  have 
produced  flow  reversal  regions  which  originated  near  the  airfoil 
trailing  edge  and  moved  upstream.  However,  this  experiment  produced 
flow  reversal  regions  which  initially  appeared  near  the  airfoil 


leading  edge  and  spread  downstream. 

It  is  shown  that  the  influence  of  unsteady  flow  parameters 
on  the  onset  of  flow  reversal  in  the  two  types  of  experiments  are 
similar.  Using  a fixed  airfoil  experiment,  the  flow  reversal 
measurements  above  the  surface  of  the  airfoil  can  be  more  easily 
and  more  accurately  obtained. 

Preliminary  measurements  of  flow  reversal  in  steady  flow 
indicated  a small  reversal  region  near  the  airfoil  leading  edge, 
which  represented  part  of  the  laminar  separation  bubble.  These 
measut ements  also  determined  the  steady-flow  stall  angle  of  attack, 
which  provided  a division  between  the  different  types  of  reversal 
phenomena  observed  in  unsteady  flow. 
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At  both  low  and  high  angles  of  attack,  the  unsteady  effects 
were  seen  to  be  caused  by  unsteady  temporal  and  spatial  quantities. 
Variations  in  the  reduced  frequency  and  the  amplitude  ratio  of  the 
oscillating  free  stream  velocity  were  observed  to  produce  changes  in 
the  relative  influence  of  the  temporal  and  spatial  unsteady 
gradients.  The  influence  of  variations  of  oscillation  frequency 
and  amplitude  were  observed  to  be  strongly  coupled,  because  of  the 
nature  of  the  unsteady  boundary  layer  equations.  The  amplitude 
effects  were  seen  to  depend  strongly  on  whether  the  amplitude  ratio, 
N,  was  near  unity.  These  amplitude  effects  had  not  been  previously 
measured  experimentally  for  the  reversed  flow  region  above  the 
airfoil  surface. 

For  angles  of  attack  below  the  steady  stall  angle,  the  flow 
reversal  response  was  close  to  quasi -steady  at  low  amplitude  ratios. 
An  increased  amplitude  ratio  greatly  increased  the  size  of  the 
reversed  flow  region  and  caused  it  to  lag  the  quasi -steady  phase. 

No  significant  unsteady  effects  we re  noted  for  variations  of 
reduced  frequency  at  this  angle  of  attack. 

For  angles  of  attack  greater  than  the  steady  flow  stall 
angle,  the  unsteady  spatial  and  temporal  gradients  cause  departures 
from  the  quasi -steady  response.  The  phase  and  spatial  extent  of 
the  reversed  flow  region  were  seen  to  be  strong  functions  of 
reduced  frequency,  amplitude  ratio,  and  angle  of  attack,  with  the 
response  normally  lagging  the  quasi-steady  flow,  except  at  high 
angles  of  attack.  The  influence  of  higher  harmonic  content  in  the 
power  spectra  of  the  free  stream  velocity  oscillations  is  observed 
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to  contribute  to  a response  which  leads  the  phase  of  the  quasi  - 
steady  flow. 

The  shape  of  the  reversed  flow  region  noted  at  various  times 

in  the  oscillation  period  lends  credence  to  the  reversal  and 

5 

separation  model  of  Sears  and  Telionis.  For  certain  combinations 
of  the  unsteady  flow  parameters,  the  observed  thickening  of  the 
reversed  flow  region  near  the  trailing  edge  may  indicate  flow 
separation.  Increased  oscillation  frequencies  appear  to  inhibit 
separation  near  the  trailing  edge.  Without  additional  data  on  the 
velocity  profiles  and  displacement  thicknesses  throughout  the 
boundary  layer,  no  further  validation  of  this  model  can  be  made. 

Although  the  split-film  sensor  proved  valuable  for  measuring 
the  onset  of  flow  reversal,  it  was  unable  to  provide  accurate 
measurements  of  the  magnitude  of  the  reversed  flow  velocity.  It  is 
recommended  that  the  velocity  magnitude  and  direction  in  the 
reversed  flow  region  be  measured  using  a split-film  sensor  with  the 
plane  of  the  split  parallel  to  the  longitudinal  flow  axis.  Also, 
another  measurement  device,  such  as  a laser  doppler  velocimeter, 
might  be  used. 

The  higher  harmonic  content  of  the  free  stream  velocity 
oscillations  was  reduced  by  adjustment  of  the  bypass  flow  of  the 
rotating  vane  device.  The  best  operating  conditions,  with  a high 
reduced  frequency  and  a small  amplitude  ratio,  provided  a second 
harmonic  rejection  of  approximately  48  db.  Other  combinations  of 
flow  parameters  used  in  the  experiment  produced  second  harmonic 
rejections  of  at  least  10  db.  It  is  recommended  that  a feedback 
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control  system  for  the  motor  drive  similar  to  that  used  by  Simpson 
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be  used  to  greatly  improve  the  harmonic  rejection  ability  of  this 
unsteady  flow  device. 
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1.  PROGRAM  1 FOR  TI-58  PROGRAMMABLE  CALCULATOR 
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Description: 

This  program  computes  the  position  of  the  probe  with  respect  to  the 
airfoil  trailing  edge  at  a given  angle  of  attack. 


Instructions: 


1.  Enter  angle  of  attack  (a) 
in  degrees 

2.  Enter  XQ,  position  of 
trailing  edge,  in  inches 

3,  Enter  Y0,  position  of 
trailing  edge,  in  inches 

4,  Enter  X-traverse 


5,  Enter  Y-traverse 


Enter 

Press 

Dis£lax 

a 

STO  00 

a 

Xo 

STO  01 

Xo 

Yo 

STO  02 

Yo 

Xtrav 

A 

X from  trailing 
edge 

B 

X from  leading 
edge 

Ytrav 

C 

Y from  trailing 
edge 

Program  Steps: 
Location 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


Code 

Bx. 

76 

2nd-LBL 

11 

A 

75 

r- 

43 

RCL 

01 

01 

95 

= 

94 

+/- 

42 

STO 

03 

03 

91 

R/S 

76 

2nd-LBL 

12 

B 

43 

RCL 

03 

03 

55 

V 

53 

( 

43 

RCL 

00 

00 

39 

2nd-C0S 

54 

) 

55 

01 

i 

00 

0 

95 

= 

75 

- 

01 

l 

95 

= 

94 

+/- 

42 

STO 

04 

04 

91 

R/S 

/6 

2nd-LBL 

13 

C 

75 

- 

43 

RCL 

02 

02 

95 

= 

42 

STO 

05 

05 

91 

R/S 

81 

RST 
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2.  PROGRAM  2 FOR  TI-58  PROGRAMMABLE  CALCULATOR 
Description: 

This  program  computes  the  height  cf  the  probe  above  the  airfoil 
surface  at  a given  probe  position  from  the  trailing  edge  and  at  a 
given  angle  of  attack.  The  airfoil  surface  coordinates  of  the  NACA 
0012  profile  are  stored  in  calculator  memories,  The  input  probe 
positions  were  calculated  using  Program  1. 


Instructions: 


Enter 

Press 

Display 

1. 

Enter  angle  of  attack  (a) 

a 

STO  00 

a 

2. 

Enter  iteration  counter 

37 

STO  06 

37 

3. 

Enter  Xy^  - from  Program  1 

X 
— 1 
m 

A 

4. 

Enter  - from  Program  1 

XLE 

B 

5. 

Enter  Yy^  - from  Program  1 

yte 

C 

AY 

Store  airfoil  surface  coordinates  below  in  memories  indicated: 
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Program  Steps: 


: at ion 

Code 

Kex 

Location 

Code 

tel 

0 

76 

2nd-LBL 

50 

75 

1 

11 

A 

51 

73 

RCL  2nd-IND 

2 

42 

STO 

52 

07 

7 

3 

03 

03 

53 

95 

= 

4 

76 

2nd-LBL 

54 

55 

5 

12 

B 

55 

53 

( 

6 

42 

STO 

56 

73 

RCL  2nd-IND 

7 

04 

04 

57 

06 

6 

8 

76 

2nd-LBL 

58 

75 

- 

9 

13 

C 

59 

73 

RCL  2nd-IND 

10 

42 

STO 

60 

07 

7 

11 

05 

05 

61 

54 

) 

12 

97 

2nd-DSZ 

62 

95 

13 

06 

6 

63 

65 

X 

14 

40 

2nd-IND 

64 

53 

( 

15 

06 

6 

65 

73 

RCL  2nd-IND 

16 

32 

X 

VA 

<■+ 

66 

08 

8 

17 

43 

RCL 

67 

75 

- 

18 

04 

04 

68 

73 

RCL  2nd-IND 

19 

22 

INV 

69 

09 

9 

20 

77 

2nd-x  > t 

70 

54 

) 

21 

15 

E 

71 

85 

+ 

22 

91 

R/S 

72 

73 

RCL  2nd-IND 

23 

76 

2nd-LBL 

73 

09 

9 

24 

15 

E 

74 

95 

= 

25 

43 

RCL 

75 

42 

STO 

26 

06 

06 

76 

10 

10 

27 

75 

- 

77 

73 

RCL  2nd-IND 

28 

01 

1 

78 

08 

8 

29 

95 

= 

79 

75 

- 

30 

42 

STO 

80 

73 

RCL  2nd-IND 

31 

07 

07 

81 

09 

9 

32 

73 

RCL  2nd-IND 

82 

95 

= 

33 

06 

6 

83 

55 

V 

34 

85 

+ 

84 

53 

( 

35 

01 

1 

85 

73 

RCL  2nd-IND 

36 

08 

8 

86 

06 

6 

37 

95 

= 

87 

75 

- 

38 

42 

STO 

88 

73 

RCL  2nd-lND 

39 

08 

08 

89 

07 

7 

40 

73 

RCL  2nd-IND 

90 

54 

) 

41 

07 

7 

91 

95 

42 

85 

+ 

92 

55 

■f 

43 

01 

1 

93 

01 

1 

44 

08 

8 

94 

00 

0 

45 

95 

= 

95 

95 

= 

46 

42 

STO 

96 

22 

INV 

47 

09 

09 

97 

30 

2nd-TAN 

48 

43 

RCL 

98 

42 

STO 

49 

04 

04 

99 

11 

11 

121 


Location  Code 


100 

09 

9 

101 

00 

0 

102 

75 

- 

103 

43 

RCL 

104 

11 

11 

105 

95 

= 

106 

38 

2nd-SIN 

107 

42 

STO 

108 

12 

12 

109 

09 

9 

no 

00 

0 

111 

75 

- 

112 

43 

RCL 

113 

00 

00 

114 

85 

+ 

115 

43 

RCL 

116 

11 

11 

117 

9- 

= 

118 

38 

2nd-SIN 

119 

35 

1/x 

120 

65 

X 

121 

43 

RCL 

122 

12 

12 

123 

95 

= 

124 

65 

X 

125 

43 

RCL 

126 

10 

10 

127 

95 

= 

128 

65 

+ 

129 

53 

( 

130 

43 

RCL 

131 

03 

03 

132 

65 

X 

133 

53 

( 

134 

43 

RCL 

135 

00 

00 

136 

30 

2nd-TAN 

137 

54 

) 

138 

54 

) 

139 

75 

140 

43 

RCL 

141 

05 

05 

142 

95 

= 

143 

94 

+/- 

144 

91 

R/S 

145 

81 

RST 

t 
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UNSTEADY  FLOW  ANALYSIS 

This  simplified  analysis  of  the  unsteady  wind  tunnel  flow 

with  identical  main-flow  and  bypass  rotating  vanes  was  adapted  from 

30 

the  analysis  of  McMichael.  It  shows  the  adjustment  required  for 
the  bypass  port  to  obtain,  ideally,  only  the  add  harmonic  components 
of  the  oscillating  velocity.  An  extended  analysis  for  the  case 
where  the  main-flow  and  bypass  vanes  are  not  the  same  size  shows 
that,  while  the  equations  are  a great  deal  more  complex,  optimiza- 
tion is  still  possible. 

The  momentum  equation  for  the  wind  tunnel  flow  may  be  written 


1“  + u ^ 
at  ax 


(f)-  y(x)“: 


where  y[x)  represents  losses  due  to  the  friction  of  the  walls.  The 


mass  balance  is: 


Q(x,t)  = A(x)  u (x , t ) 


where  A(x)  is  the  local  cross-sectional  area.  Substituting, 

1 9Q  + Qi  L ( 1 \+  QMx)  = _ J_  (l \ 

A(x)  at  2 ax  \Az(x)/  A2 (x)  ax  \ p / ' 

The  wind  tunnel  stations  and  mass  flows  are  defined  in  Figure  7. 
Integrating  between  arbitrary  stations  m and  n yields: 

<PnrPn>  = “„,n  S + S*  1)2 


where: 


“mn  = p/ 


ATxT 


Kmn = \ Miry], + p/  dx- 


r ~ ^ 
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Knowing  that  P = P = Pa  (atmospheric  pressure)  and  assuming  quasi- 
steady,  high  Reynolds  number  flow,  the  pressure  differences  at  each 
station  are: 

P.-Pi  ■ % Q,  + I,2 

P.-P.  ' % «,  + V 

P,-P«  ‘ M,  + Q2>  + ks6  (Q,  +Q2)2 

p,-p.  ■ <vq2> + tvu1 

al  so: 

P - P = C Q 2 

2 3 23  W1 

P,-PS  = C„5  V 
p,-ps  = Cas  Q** 

where  Cmn  is  the  pressure  drop  coefficient  across  the  indicated 
mechanical  devices.  Also: 

P -P  = -P r (fan  pressure  rise). 

6 7 T 

Then,  McMichael  writes: 

(P,-Ps)  - <Pa-Ps)  = 


Now  define 


a.  = (a,,  + a ) 

1 '12  3i)  ' 

K,  = (K  + K +C  ) 

1 12  34  23 


Fi  = C,s 


F2  " Ca5* 


Combining  these  terms  gives 
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F2Q22  = olQ1  + (K!  + Fj )Qj2* 


The  adjustable  bypass  ports  can  be  opened  a distance  L and  essen- 
tially add  a resistance  K2(l)  so  that 

Pa-P5  = (K2  + F2)Q22  . 

Thus 

(Pa-P5)  = (K2  + F2)Q22  = ctjQj  + (Ki  + Fx  )Q| 2 . 

Note  that  ot2Q2  was  considered  negligible  because  of  the  short 
distance  between  the  bypass  port  and  bypass  vanes.  For  the  fan 
duct  exit 


(P5-Pal  = (Ps-Ps)MP6-P,)  + (P7-P,) 

" Ks  +a?a><Q1  +Q;)  + (Ks6  + KjJIQ,  + Q;);  - Pf 


define 


K3  = KSS 


+ K_ 


78 


a 


56 


78 


Collecting  the  previous  relations  yields 

'(k2-*f2)Q2s 

Pa-P5  = ■ (K,  + F.JQ,2  + 0..Q, 

■-K,(Q,+Q2)!  - a.IQ.+Q,)  ♦ P f- 


(B-l) 

(B-2) 

(B-3) 


A reference  state  for  the  tunnel  flow  may  be  defined  using  the 
time  average  of  the  vane  pressure  drop  coefficients  and  mass  flows 
at  a given  mean  tunnel  velocity.  Then, 


where 


also, 


F,<»>  ' F„  * f,(t) 
F,(*>  ' F20  + f2(t) 
Q,(t)  - Q10  * q,(t) 
Q, ( t)  = Q,.  + q,(t) 


F7W=  F, 
F7W-  F, 

57ftT  * Q, 

Q7!tT  - Q. 


f,U)  ■ f2(t)  = q,(t)  = q2(t)  = o 


P-  = Pf  + AP* 
i 1 n < 


-Hi-  ^ _3l_  ^ £ « 1 

Q ^ Q ^ Pf  <<  1 • 
W10  V20  To 


The  reference  state  (time  average)  of  equation  s (1),  (2),  and  (3) 


<VF,.>Q2.!  - ^Ki  + Fio  >^io2 

<K,tF,.M„,2  ' Ff  “ f3(Q10  + Q20  )z  • 


(EM) 

(B-6) 


Equations  (B-l),  (B-2),  and  (B-3)  are  then  linearized  about  the 
reference  state  (second-order-terms  dropped): 

F(F2  + F20  )Q20Q2  + f2Q2(Jz  = 2(K1*F1,)Qltq,  + f1Q112+a]q|  (B-6) 
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2(K2  + F20)Q20q2  + f2Q;02  = APf  - a3(q1  + q2)  - K$(Q10  + Q20  )(q4  + q2) . 

(B-7) 

For  the  simplified  case  where  the  main-flow  and  bypass  vanes  are 
identical,  equations  (B-4)  and  (B-5)  show  that  F = F 5 F . To 
balance  the  reference  state  flow  rates,  K2  is  adjusted  so  that 
Kj  = K2,  so  that  Q10  = Qq/2 . Then  define: 

Ri  = Ki  + Fio 

Ra  = K2  + F20 
r2  = Rj/R2 

Q2o  = rQio- 

The  equations  (B-6)  and  (B-7)  become: 

2rR2Q10q2  + f2r2Qio2  = 2R1Q10q1  + f ^ + c^q,  (B-8) 

= APf -a3^i +^2)  -2K3Q10O+r)(qi  + q2) 

for  simplification,  APf  may  be  approximated  as 

APf  * B(q,  +q2) 

where  B is  a constant.  Equation  (B-8)  may  be  non-dimensional ized 
using 

Rj  = R2  = R' 
f,  = f ,/2R' 

f2  = f 2 /2R ' 

5,  - q,/(Q0/2) 


q,  ■ q,/(Q0/2) 


which  gives 


2q2  + 2f2  - 2q,  + 2f , + 


23  / ^ A \ 2cy  q / a a 

- RQ-  (q!  - q2)  - rq-  (qt  - q2) 


(5i  + q2) - (B-9) 


Now  define  the  time  constants 


and  also  let 


RQ0  ’ 


r _ _G_ 
^ nrv 


Substituting  gives 


T!  t3  q,  + [Tj  (1  -G)  + 2x3]^  + (l-2G)q1 

• • 

A A A A 

■ G<Vf2>  - <,(f.  -M- 


The  damping  ratio  for  this  second  order  system  is 


^ = (l-G)xi  + 2t3 


2 A(1  -G~)x ! x , 


and  the  resonant  frequency  is 


R(L  /1-2G 


( B— 1 0) 


. /T2G  _ 

O r j 


It  can  also  be  shown  that  these  small  vanes  of  set  A ideally 
generate  only  even  harmonics  in  the  velocity  field.  Equation  (B-10) 

/\  A 

shows  that  the  forcing  function  is  (f3  -f2).  Figure  B-l , also  from 
McMichael , shows  the  angular  dependence  between  the  two  sets  of 
vanes  to  be 
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2<J>j  = a)t  + tt/2 


2<}>2  = u)t  - ir/2. 


For  a given  vane  resistance  function 


f(<j>)  = -£  + X)  an  cos  (2n<j>) 
* n=l 


it  can  be  seen  that 


+ ^2^2)  " + S 


? ■ t_j  a cos  (n(wt  + ir/2)) 

c n=0,2,4, . . . n x ' 


and 


f ,(<!>,)  - ?2 (4>2 ) = 2 an  cos  (n(wt  - v/2}) 

n-1  J 3 ) 5 ) a a a 

A A 

The  forcing  function  (fx  - f z)  in  the  dynamic  equation  thus  ideally 

only  contains  odd  harmonics.  This  was  difficult  to  determine 

experimentally  since,  for  the  small  vanes,  usually  only  the  second 

harmonic  was  of  a large  enough  amplitude  to  be  detected  by  the 

Fourier  analyzer  on  a logrithmic-magnitude  power  spectrum  (maximum 

scale  of  80  db).  The  equations  cannot  be  so  easily  simplified  for 

the  situation  where  the  main-flow  vanes  and  bypass  vanes  are 

different,  because  F f F and  K 1 K . The  first  line  of 

10  20  12 

equation  (B-8)  may  then  be  written  (including  the  assumption  for 
APf): 

q2  = A'f,  - B ' f 2 + C'q,  + D'q,  (B-ll) 

where 
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+ 2i^[2K>Q“.(1+r)r’-s,"!]  (B-12) 


The  effect  of  the  variable  bypass,  through  R2,  is  seen  to  affect 

i 

each  term  of  equation  (B-12). 


/ 


